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Abstract 

The phenomenon of super-radiance (Dicke effect, coherent spontaneous radiation 
by a gas of atoms coupled through the common radiation field) is well known in quan- 
tum optics. The review discusses similar physics that emerges in open and marginally 
stable quantum many-body systems. In the presence of open decay channels, the in- 
trinsic states are coupled through the continuum. At sufficiently strong continuum 
coupling, the spectrum of resonances undergoes the restructuring with segregation of 
very broad super-radiant states and trapping of remaining long-lived compound states. 
The appropriate formalism describing this phenomenon is based on the Feshbach pro- 
jection method and effective non-Hermitian Hamiltonian. A broader generalization 
is related to the idea of doorway states connecting quantum states of different struc- 
ture. The method is explained in detail and the examples of applications are given 
to nuclear, atomic and particle physics. The interrelation of the collective dynamics 
through continuum and possible intrinsic many-body chaos is studied, including uni- 
versal mesoscopic conductance fiuctuations. The theory serves as a natural framework 
for general description of a quantum signal transmission through an open mesoscopic 
system. 
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1 Introduction 



More than half a century ago Robert Dicke published an article where he pre- 
dicted the phenomenon of super-radiance [1]. This phenomenon was indeed 
experimentally observed in an optically pumped gas [2] and later studied in de- 
tail and widely used in quantum optics [31 H] as well in physics of plasmas and 
beam physics [S]. The super-radiance was also observed with the free-electron 
laser [5] and with atomic matter waves, including Bose-Einstein condensate [7]. 
The quantum-mechanical essence of super-radiance in the original formulation 
is relatively simple. Two-level systems ("atoms" or "qubits") can be coupled 
even in the absence of their direct interaction. The coupling proceeds through 
the common radiation field. The migrating resonant excitation combines the 
atoms in a collective capable of releasing the burst of radiation as a whole. The 
phenomenon therefore is called coherent spontaneous radiation that sounds as 
an oxymoron since traditionally the spontaneous radiation was considered as 
being incoherently produced by individual atoms. 

In the end of the nineteen eighties it was understood that the physics un- 
derlying super-radiance is much more general and can find broad applications 
in various regions of the quantum world. The necessary ingredient is the set 
of unstable quantum states with the same exact quantum numbers. The states 
should be capable of decaying into a common decay channel (or few channels). 
This creates an analog of their coupling through the common "field of radia- 
tion" , or common continuum; the coupling becomes effective when the widths 
of individual adjacent states grow and overlap. The system of overlapping res- 
onances undergoes a sharp restructuring forming super-radiant states, one for 
each common open channel. In the ideal version of quantum optics, when the 
wavelength of radiation significantly exceeds the size of the system, the identical 
atoms have identical levels and widths which amplifies the effect. In a general 
case, the levels are not degenerate, and the super-radiant transition can happen 
only when the typical widths are of the order of typical spacings between the 
unstable levels. 

The mechanism of super-radiance is omnipresent being simply a consequence 
of quantum mechanics. It was seen in the problems of quantum chemistry, 
atomic and molecular physics. In a pure form it was discovered in nuclear 
physics in consideration of the decay of giant resonances [8 when the numerical 
simulation clearly indicated the appearance of very broad resonances for each 
open decay channel when the continuum coupling exceeds a critical strength. 
At the same time the remaining resonances are getting more narrow signaling 
the existence of trapped states with long lifetime. 

The physical and mathematical analogy of this behavior of a set of reso- 
nances to the super-radiance phenomenon was first indicated in [5J [TU] . The 
starting point there was quantum chaos, namely a question of the influence of 
the continuum coupling on the statistical distribution of the level spacings. It 
turns out that the super-radiant behavior survives the intrinsic dynamical chaos 
being mainly determined by the character of decay. The spectral statistics, typ- 
ical for canonical ensembles of random matrices, evolve into the statistics of 
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complex energies which are at the same time the poles of the scattering matrix 
in the space of open channels. The statistical properties, in turn, are instru- 
mental in defining the fluctuations and correlations of cross sections of various 
reactions proceeding through the channels open in the energy interval under 
study. The phenomenon of super-radiance reveals itself in a specific two-cloud 
distribution of energies and widths of resonances in the complex plane [11] [T^] ■ 
The adequate mathematical formalism for the description of this physics is 
provided by the efFective non-Hermitian Hamiltonian |13j suggested in nuclear 
theory. The method goes back to Feshbach [HI [15] and Fano [16] extending 
and improving (eliminating not well defined channel radii) the boundary con- 
dition method originated by Kapur and Peierls [T7] , Wigner and Eisenbud [TS] , 
and Teichmann and Wigner |19| . Similar ideas can be found already in much 
earlier notebooks by Majorana [20]. The classical review articles are available 
concerning this period of development [H] |22l [23]. The equivalent method 
was independently developed and applied, among other problems, to quantum 
scattering theory by the outstanding Russian-Israeli mathematician M.S. Livsic 

This approach was described in the well known book [26 although the fact 
that the super-radiant physics directly follows from this Hamiltonian was un- 
derstood only later. The efFective Hamiltonian, in general energy-dependent, 
can be derived by a standard projection method when the part of the dynamics 
corresponding to the continuum channels is formally eliminated and the whole 
Schrodinger equation is projected onto the internal dynamics. In this part of 
Hilbert space the non-Hermiticity appears as a result of the presence of open 
channels so that the eigenstates above threshold become unstable resonances, 
and the eigenvalues of the Hamiltonian are complex energies £ = E — (j/2)r, 
where F is the width of the resonance. The anti-Hermitian part of the Hamil- 
tonian has a specific factorized structure as follows from the unitarity of the 
dynamics in the full space [27] . In the regime of strong continuum coupling this 
part plays the main role and the rank of the corresponding matrix, equal to the 
number of open channels, determines the number of super-radiant states. The 
Hermitian part of the effective Hamiltonian is responsible for energy shifts of 
resonances and bound states; by its nature it is similar to the Lamb shift but 
in this case of collective origin. It is useful to stress that the non-Hermitian 
character of the Hamiltonian here is not related to any approximation as long 
as all channels are accounted for. 

In our review, the formalism of the effective Hamiltonian will serve as a 
foundation of one of the most effective approaches to the description of open 
many-body systems. This is especially important for modern nuclear physics 
where the center of interest moved to nuclei far from stability. Such nuclei are 
often only marginally stable being close to the drip lines which delineate the 
limits of existence for nuclear systems. The proximity of the continuum — both 
actual decay channels and virtual intermediate states — determines the specific 
physics of such systems. The continuum shell model [551 IMl [33 [SI] developed 
in the framework of the effective Hamiltonian allows one to describe in a unified 
approach bound states, resonances and reaction cross sections. In distinction to 
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other perspective approaches to loosely bound nuclei, such as the Gamow shell 
model, that uses the analytical continuation to the complex plane [32 123 El] , 
here one can consider all processes at real energy of the reaction. The energy 
dependence of the effective Hamiltonian reflects, in particular, the existence of 
channel thresholds, and the proper behavior of the reaction amplitudes at the 
thresholds is a noticeable and practically important feature of the continuum 
shell model. 

Nuclear physics is only one, even if essential, area of applications of the ideas 
of the effective Hamiltonian and super-radiance. In fact, this is a convenient way 
to describe dynamics of any quantum system in the energy region including both 
bound states and continuum. Atomic autoionizing states [3S] and intermolecular 
processes |36| were long ago shown to reveal a similar behavior. Furthermore, 
we can expect the same trends to manifest themselves in physics of intermediate 
and high energy [37l |38l |39] . The role of open channels in various situations 
can also be played by the doorway states [10] which are forming a bottleneck 
connecting the states of different nature [H]; this is another generalization of 
the ideas of super-radiance. 

Important applications of the super-radiance and related theory can be found 
in mesoscopic physics. The universal mesoscopic conductance fluctuations [42) 
share significant similarity with Ericson fluctuations [121 [H] typical for nuclear 
reactions observed in the energy region of overlapping resonances [45] . The 
detailed analysis of correlations and fluctuations of cross sections for different 
channels can shed new light on the fundamental problem of the signal transmis- 
sion through a quantum system, from simplest qubits to very complex molecules 
and artificial structures with interfering propagation paths. The first results 
modifying the Ericson theory by the effects of super-radiance and dependence 
of certain observables on the degree of internal chaoticity appeared recently in 
the literature pS', '57! , along with the applications to ordered and disordered 
nanostructures [49 1 iSOj . Vice versa, the intrinsic chaos is strongly modified by 
the openness of the system. 

Our review is devoted to the manifestations of the super-radiant mechanism 
outside of quantum optics. The presentation is organized as follows. We start 
with illustrations of the phenomenon of super-radiance in simple cases. Chapter 
2. In Chapter 3 we introduce the projection formalism, derive the effective 
Hamiltonian and show how the generalized super-radiance emerges as a result of 
the strong coupling of the intrinsic states to the continuum channels. As the first 
application of the effective Hamiltonian we consider the nuclear continuum shell 
model. Chapter 4 contains examples showing this approach at work. In Chapter 
5 we show the extension of the approach to doorway states and such collective 
excitations as giant resonances. Other examples of the same approach. Chapter 
6, will illustrate applications to different areas of physics. We devote Chapter 7 
to the problems of quantum chaos in open systems and their manifestations in 
spectral statistics and reaction cross sections. Finally, Chapter 8 will discuss the 
rapidly developing area of quantum signal transmission from the same general 
viewpoint. The conclusion will be used mainly to list the interesting problems 
which are still waiting for their explorer. We will try to explain the main physical 
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ideas along with their mathematical formulations but we will not be able to go 
through detailed calculations referring instead to the original literature. 
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2 Origin of super-radiance 



2.1 Original idea 

For a simple start we will go back to the original idea of optical super-radiance 
[T] and stress the elements which later will be generalized to similar phenomena 
in various subfields of physics. Let us consider an idealized system of N iden- 
tical two-level atoms confined to the volume of a size much smaller than the 
wavelength of radiation that corresponds to the transition between the levels of 
each qubit. We can describe the intrinsic states of such a system with Na 
atoms in the ground state and Ni, ^ N ~ Na atoms in the excited state with the 
aid of corresponding creation, and b\ and annihilation, a and b, operators as 

\^iNa,N,)) ^ \Na,N,) = {bY' (a^riO), (1) 

where for simplicity we assume that the atoms are bosons and |0) denotes the 
vacuum state. 

Neglecting the radiation from the excited level, we have stable intrinsic space 
spanned by the states Q. It is useful to apply the Schwinger representation of 
angular momentum |51| with the aid of individual spins Si — 1/2 corresponding 
to each qubit and the total spin S = ^ - of the entire system that is usually 
called here the Bloch vector. The identification of quantum numbers is obvious: 

\^>{Na, Nt)) ^ \SM), ^{Na + Nt)^ ^N, M ^ ^ ^ {N, - iV,), 

(2) 

or, in terms of original atomic operators, 

\SM) = [{S + A/)!(5 - (5t)S+M(„t)S-M|Q^_ 

The Hamiltonian of the closed system is trivial, 

H = uJJ2^^z=^M, (4) 

i 

where the origin of the energy scale is set at the middle between the excited and 
ground levels of an atom. 

Now we make the excited atoms unstable introducing the radiation field. If 
A* and A are, in appropriate units, the amplitudes of radiation and absorption 
of the photon, respectively, this part of the Hamiltonian is 

H' = A*a^b + Ab^a = A*S-+AS+, (5) 

where, using the long wave limit, we neglect the phase differences between the 
atoms. Each act of radiation brings M by one step lower. If all atoms are 
initially excited, S — M = N/2, the Bloch vector is at the north pole of the 
Bloch sphere, and only spontaneous radiation of each atom is possible with 
the total probability \A*{S,S - 1\S^\S,S)\^ = \A\^2S = \A\^N. This is an 
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incoherent process of independent radiation by excited atoms with the intensity 
proportional to their number. 

However, the same Hamiltonian ([5| describes also the processes of interac- 
tion between the atoms though virtual radiation and absorption of quanta. In 
this idealized model, such processes are amplified by the exact degeneracy of 
atomic levels all tuned in resonance. While the spontaneous radiation brings the 
average level of excitation down and the Bloch vector deviates from the vertical 
position, the amplitude of radiation grows through the virtual coupling between 
the radiators that becomes possible due to the devastation of large-M states and 
opening possibility of the virtual absorption. The maximum is reached at M = 
(the Bloch vector in the equatorial plane) being proportional to |ApS'(S' + 1). 
Here we have a coherent spontaneous radiation — super-radiance — with the 
probability oc iV^. Such a burst of spontaneous radiation was observed in opti- 
cally pumped HF [2] and later studied in numerous experiments. The underlying 
dynamics of the Bloch vector can be presented also in a language of classical 
mechanics [5^ . 

Of course, the practical implementation is not identical to this oversimplified 
model. The different spatial positions of the atoms bring in the phase differences 
between interfering radiators. (In fact, the seemingly opposite limit of the size 
large compared to the wave length allows one to measure the so-called collective 
Lamb shift of the radiation line [53] ; even a very small shift of this nature can 
limit accuracy of atomic clocks [M]). The local nonuniformity and collisional 
broadening lift the exact degeneracy of atomic levels. Direct, for example dipole- 
dipole, interactions make the system different from the collection of independent 
radiators, the extension of the model still allows an exact solution [55l [56] . 
However, even if not in such a pure form, the super-radiance appears as an 
important physical effect with abundant applications; the internal states of the 
radiators are effectively coupled through the common radiation field. 

For our goal, the main lesson is that the quantum states, along with possible 
direct mixing by interactions, can be coupled also through the continuum of open 
decay channels. Since the continuum coupling determines the width F, or the 
lifetime r ~ ?i/F, the states become quasistationary and can be characterized 
by a complex energy, 

£ = E-'-T. (6) 

Similarly to standard perturbation theory, the efficiency of coupling is deter- 
mined by the ratio of the coupling strength to the energy spacing between the 
coupled states. If the width F is of the order of, or exceeds, the spacing between 
real energies E, coupling through the continuum turns out to be effectively 
strong. This transition to the regime of overlapping resonances supplies the 
system by the new collectivity - the possibility of coherent decay. 

2.2 Transport through a lattice 

The idea of generalized super-radiance works for many open quantum systems. 
We illustrate this by a seemingly very different example [13_ taken from con- 
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densed matter physics, where we will see the similar character of emerging 
effects. Let us consider the simplest one-dimensional crystal with identical cells 
numbered 1, 2, N. A particle has levels at energy e in all cells and the hopping 
amplitude v to adjacent cells. The edge barriers allow the coupling to environ- 
ment so that the corresponding decay widths are 7^ and 7^^. The system can 
be described by an effective non-Hermitian Hamiltonian with matrix elements 
in the site representation 

which describe the decay to outside with the help of complex energies Q at 
the edge sites. In the absence of decay, the internal levels would be split by the 
hopping and form a crystal band of Bloch waves labeled by the discrete wave 
vectors q — 1,...,N. The band is characterized by corresponding energies (we 
set e = 0) 

eq^2vcostpq, ipq= j^^^ Q, (8) 
and the eigenvectors |g) with components in the site representation 



{n\q) = Y sin(n^g). (9) 

In this basis, the Hamiltonian ([7]) is transformed to 

with the factorized matrix elements of the anti-Hermitian part, 

Wqq' = Sin^q Sm^q,[jL + + Jr]- (H) 

As the eigenstates are delocalized, each of them acquires a width; the maxi- 
mum width appears in the middle of the energy band. The intrinsic hopping 
interaction is now substituted with the anti-Hermitian interaction through the 
continuum. The quasistationary states can be found by diagonalizing the effec- 
tive Hamiltonian (11). The complex energies ^ are given by the roots of the 
secular equation 

l + lijL+lR)P+{£)+lLlR[Pl{£)-Pl{£)]^0, (12) 



where 



(13) 



In all cases the total width of all resonances J^q to be equal to the 

trace, 7^+7^', of W. Two possible limiting situations clearly come out. At weak 
continuum coupling, the individual Bloch states give rise to isolated resonances 
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which can still be labeled by the wave vector q. Neglecting the last term in eq. 



^a^^,-^^,^ r, = W^,„ (14) 



(12), we find the widths as the diagonal elements of W , 

i _ 
2' 

which, by use of eq. (|9]), can be presented in an obvious form 

r,-7L|(i|g)P+7fil(iV|g|)P- (15) 

The decay here serves as an analyzer singling out the components of the qua- 
sistationary state matched to the specific decay channel, left or right in this 
model. Such arguments are usually used as a justification [57' for identifying 
the width distribution in a compound nucleus with the distribution of the sin- 
gle components of the compound nuclear wave function. For reduced neutron 
widths of low-lying neutron resonances with Gaussian distribution of basis com- 
ponents this gives rise to the Porter-Thomas distribution (PTD) [SS]. However, 
this result and underlying logic are valid only in the lowest approximation with 
respect to weak continuum coupling of non-overlapping unstable states. 

In the opposite limiting case of strong continuum coupling, the lifetime with 
respect to the irreversible decay outside is shorter than for the hopping between 
the sites, Jl,r ^ v- In this limit (and N :S> 1), one can neglect the band 



splitting Eq in the denominators of P± in eq. (13), so that P_ vanishes while 
P+ l/{2£). Then two roots are shifted along the imaginary energy axis 
defining 

£l,r = - 2 (16) 

This limit of short-lived collective states is an analog of optical Dicke super- 
radiance. (If 7l and significantly differ, the system will reveal two super- 
radiant transitions j50|). In distinction to usual collective states in a many-body 
system, here the collectivization proceeds through the widths of resonances, all 
Bloch waves radiate coherently being coupled through the decay amplitudes. 
In the g-space, the super-radiant states are in the middle of the spectrum, 
while in the site representation it is trivial to find that they are located at the 
edges. Therefore, at strong continuum coupling it might be useful to go to the 
doorway representation [13[ that will be discussed in the corresponding section 
from a more general viewpoint. In the one-channel case, the transformation 
of the S'-matrix to such a doorway form can be made explicitly as was shown 



long ago [59]. In the limit (16), remaining N — 2 intrinsic states are darJc 
being decoupled from the continuum. It is easy to find the correction to eq. 
(16), where the long-lived (trapped) states acquire small widths, cx {v/jl,r)'^, 



through the intermediate coupling to the doorway states directly coupled to 
the continuum. We can also mention that the disorder in the closed system 
may even amplify the effects of the super-radiance [55] through the formation 
of localized states more easily decaying outside. We return to the more detailed 
discussion of quantum transport in Sec. 8. 

The common features of our two examples - optical super-radiance and sig- 
nal transmission through a chain of radiators - are generic for all problems with 
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open quantum systems. Wc have here a kind of phase transition from inde- 
pendent unstable states to collective radiation (or any process developing with, 
or through, open continuum channels). A typical control parameter regulating 
this transformation is the relative strength of continuum coupling with respect 
to the level spacing in the closed system. Below we introduce the universal 
formalism adjusted for the consistent description of such situations. 
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3 Open quantum system 



3.1 Projection formalism 

Following [T3|[T3[23 we consider a general quantum system that has both bound 
and continuum states. The states |c; E) in the continuum will be labelled by 
total energy E and the quantum numbers of the channels. The channel, c, is 
characterized by the kind(s) of the fragments in the continuum, their intrinsic 
states, and all discrete quantum numbers (spins, isospins, parity...). Every 
channel c opens at its threshold energy E'^. The whole dynamics in the Hilbert 
space of the system is governed by the Hermitian Hamiltonian H — , while 
the reactions a — & between various channels, including the elastic one, b = a, 
can be described by the scattering matrix in the channel space, S^°'{E), which 
is unitary, — S^^. 

Using the stationary quantum mechanics we have to solve the Schrodinger 
equation for the wave function at real energy, 

H^E = E^E, (17) 

where energy E can belong to the discrete spectrum (and then has to be deter- 
mined by the solution of the equation) or to the continuum, if some channels 
are open having their thresholds below given value of E. In the coordinate rep- 
resentation, those two cases have different boundary conditions. The threshold 
positions are also to be determined self-consistently by the Schrodinger equation. 

It is always possible to formally subdivide the Hilbert space of the prob- 
lem into two (or more, see for example [60|) complementary parts using the 
expansion of the unit operator into arbitrarily selected orthogonal projectors, 

i = r + Q, rv^r, qq = q, vq = qv = o. (is) 

Then the Hamiltonian H, as any other operator, contains four classes of matrix 
elements, 

H = Hqq + Hq-p + H-pQ + H-p-p, (19) 

while the wave function has two parts, '9 = V'^ + Q'S', which satisfy the coupled 
set of equations, 

[Hqq ~E]Q^ = -HqvV^, [Hvv -E]^^ = -HpQ Q^. (20) 
Formally eliminating the part V"^ of the wave function, 

V-ii = (^E - Hpp^ ^H-pqQ-^, (21) 
we find the equation projected into the class Q, 

HiE)Q^ ^ EQ^, (22) 
where we have to deal with the effective energy-dependent Hamiltonian 

= Hqq + Hqv ^ \ HvQ. (23) 
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It is useful to stress that the formal results here and below do not require any 
explicit expressions for the projection operators although it can be done in many 
ways |15j . both in Hilbert space and in configuration space. The flexibility of 
the approach is its important advantage. Working in the Hilbert space we can 
avoid problems related to the dependence on the coordinate radius of a channel 
and explicit matching of the solutions at such a surface [M] . 



The formal solution ( 21 ) has to be complemented with rules of working with 
the singularities emerging when H-p-p has eigenvalues at energy E. Although 
the projection procedure can use any subdivision of space, we will assume that 
the part V may contain the continuous spectrum |c; E), so that the second part 



of eq. ( 20 ) has a solution of the homogeneous equation, and in any channel c 
the singularity appears at energy above the corresponding threshold {E > E'^). 
The standard rule will be applied that real energy E is to be considered as a 
limiting value from the upper half of the complex energy plane, iS^"*") = _E + iO. 
With this rule, those states of space Q which will turn out to be coupled to open 
channels in V , will acquire the outgoing waves and become unstable. Then the 



matrix elements of the effective Hamiltonian ( 23 1 can be written more in detail 
as 



■Hi2 ^Hi2 + Y^ j {dr' dT")dE'dE" 



x{l\HQ-p\c',r',E')( \ ) {c",r",E"\H-pQ\2). (24) 

Here we use an arbitrary orthonormal basis in the Q-spacc with the basis states 
labeled as |1), |2),... The differentials [dr' dr") symbolically denote kinematic 
variables in the channels c' ^c" , apart from explicitly indicated energy. 

It is always possible |64| to find the eigenchannel basis in the 'P-space that 
brings Hpp to the diagonal form with eigenvalues E'{c; r) and the Hamiltonian 



( 24 ) simplifies to 



Hi2 = Fi2 + ^ l{dr')dE'{l\HQp\c,r',E')^^^y^^^^^{c,T\E'\HpQ\2}. 

(25) 

For channels closed at energy E, there is no singularity, and the superscript (+) 
can be omitted. For open channels we use the usual identity valid for real x, 

^ P.v.—^ iTTS{E-x), (26) 



Ei+t -X ' ' E~x 

where P . v. is a symbol of the principal value. This leads us to the "canonical" 
form of the effective Hamiltonian, 

H^H- '-W. (27) 

The Hermitian part here is 

H = H + A{E), (28) 
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where A{E) is an off-shell integral in the sense of the principal value, 
Ai2(i?) = P.v.^ j {dT')dE'{l\HQr\c,T',E') ^-±^^{c,t' , E'\HrQ\2), 

(29) 

while the anti-Hermitian part is given by 

W,2iE) = 27r Al{E)Af{E), (30) 

c(opcn) 

and we introduced the effective amplitudes A'i{E) as matrix elements of -ffg-p 
taken on-shell, E' = E, with all kinematic factors coming from the integral 
included into their definition. It is important to notice that the real (dispersive) 
part A of the effective Hamiltonian, sometimes called the energy shift, includes 
the sum over all channels, open and closed, while the imaginary (absorptive) 
part W contains only contributions of the channels open at given energy. The 
kinematic factors (the density of states in the continuum) guarantee that any 
amplitude Al{E) vanishes at the threshold energy E'^. By construction, these 
amplitudes are combined into a factorized structure that is a consequence (and 
a mathematical reason) of the unitarity of the scattering matrix |57]. We have 



already seen this factorization in a specific simple model, eq. ( 11 ) 



3.2 Effective Hamiltonian and reaction theory 



The effective Hamiltonian (27 1 allows us to consider the properties of an open 
system in a very general fashion. The advantage of this formulation is in the 
discretization of the problem since we can assume that the internal part of the 



Hamiltonian, Hqq, by construction of our subdivision (19) has only a discrete 
spectrum, so that the resonances can appear only due to the continuum coupling, 
H-pQ and Hq-p. This "natural" discretization should not be confused with 
the discretization of the continuum used frequently in the numerical studies of 
various reactions [65] : here we are still working within the exact formulation of 
the theory but the whole dynamics is projected onto the "intrinsic" part of the 
total Hilbert space. 

For a system with A'^ as a (usually very large) dimension of the intrinsic 
space and M open channels, the amplitudes can be considered as rectangular 
matrices N x M which we can label A, so that the matrix of the operator W, 



eq. ( 30 1 , can be presented as a product 

VK = 27rAA^ (31) 

In a time-reversal invariant system, the amplitudes A^ can be taken as real 
numbers; then A^ = A^. This factorized structure facilitates the task of estab- 
lishing the relation of this approach to the formal theory of reactions [SBl [57] . 
On the other hand, this is a mathematical key to super-radiance. 

We can introduce the propagator at complex energy £ for the effective Hamil- 
tonian Ji as 

5{£) = (32) 
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where physical energy axis is a limit £ — ?> i?*^"*"^. The analog of the Dyson 
equation that connects Q with the propagator G for a closed system, 



G(£) = ^, (33) 

can be derived algebraically with the help of the factorized structure of W: 

Q{E) = G{E) - l G{E)A A^G{E). (34) 

2 ^ ^ l + {i/2)K{E) 

Here and later the hats mark Ad x M matrices in the channel space; in particular, 
the matrix 

k{E) ^ A^G(£)A (35) 

describes the full propagation inside the system from the entrance channel to 
the exit channel. The series of such propagations summed over in the operator 



[1 + {i/2)K] ^ of eq. ( 34 ) corresponds to the full process of the reaction through 



the system by an infinite number of intermediate excursions from the intrinsic 



space to the continuum. The matrix ( 35 ) substitutes in this theory the standard 
i?-matrix [HI [Ml [S3]- 

The part of the amplitude for the reaction b a that proceeds through the 
states of the subspace Q is now an element {TqY^ of the M x M matrix 

fQ{E)=A^g{E)A = k{E) . (36) 

l + {i/2)K{E) 

The total reaction amplitude includes also the direct processes determined by 
the dynamics in the 7^-space without entering the intrinsic state, 

f^fr+TQ, (37) 

where E — ?► E^-^K The item T-p is responsible for what is called potential scat- 
tering. The complete scattering matrix can be written as 

S{E) = 1 - iT{E) = -ifv + ^ ~ ('/^^^(^^ , (38) 

l + (^/2)i^(f)' ^ ' 

and the cross sections in open channels are defined by the corresponding matrix 
elements of the ^-matrix at physical energy, with the interference of direct 
and resonance scattering taken into account. The amplitudes A also contain 
phases of the potential scattering at the entrance and exit of the channels. 
The resonance part, closely related to the intrinsic structure of the system, will 
capture our main attention. In the absence of direct processes, this part is 
evidently unitary by itself. 



The amplitudes Tq{E), according to the first equality (36), have the poles 
in the complex plane of energy E coming from the poles of UJc), which are the 
eigenvalues ^ of the effective Hamiltonian Ti. Those eigenvalues determine the 
"compound-nucleus" resonance spectrum of the system. We use here this term 
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only to express the fact that the wave function of such a state has components in 
the Q-space. Whether a resonance indeed agrees with the notion of compound 
nucleus in the standard meaning of this concept can be decided by the com- 
parison of the resonance lifetime t ^ Ti/T with the characteristic equilibration 
time inside the system. As stated by Feshbach [TS], "the effective Hamiltonian 
method is most convenient for nuclear reactions in that it differentiates ... be- 
tween compound and direct nuclear reactions". On the real energy axis, the 
amplitudes have branch points at the channel thresholds. 



3.3 Emergence of resonances 

We start with a simple example of a resonance that emerges as a result of 



interaction between the two subspaces ( 18 ). This will illustrate a direct relation 



between the method under discussion and the textbook scattering theory. The 
formalism, apparently first used by Fano [TB] in application to the autoionizing 
states of an atom, immediately follows from the general projection method. In 
fact, this approach is similar to the consideration of the spreading width of a 
collective {^^bright") state interacting with the background of secondary states 
and in such a form it was suggested by Rice |68l . This analogy will reappear in 
the discussion of the doorway states in Sec. 6. 

If we consider only a small number of special intrinsic states, it is convenient 
to project the entire dynamics into the external 7^-space. Eliminating the part 
Q'^ of the wave function, we come to the equation 

[Hrv - E]V^ = ffpg ^ Hq-pV^. (39) 

Here energy E belongs to the continuous spectrum, and the homogeneous equa- 
tion, [H-p-p — E'jV^ = 0, has solutions ip'^, corresponding to the asymptotic 



channels c included in the P-space. The right hand side of eq. (39 1 describes 
the intrinsic states with the same conserved quantum numbers which are cou- 
pled to the continuum and, if above thresholds, can reveal autoionization. As 
a demonstration of principles, we consider here only the case of one intrinsic 
state $ with unperturbed energy Hqq = e coupled to a single channel [if not 
all open channels are included into the 7^-space, we have to consider the limit 



E — T' in the denominator of the second term in eq. (39)] 



The solution of eq. (17 1 is a superposition 



<iJE = a{E)^ + J dE'b{E,E')tjjE', (40) 

a direct analog of the spreading of a bright state in the discrete spectrum. 
Omitting from our notations the dependence of a and b on running energy E, 
we obtain the set of coupled equations for the amplitudes a and h{E') in terms 



of the continuum couphng amplitudes A{E') = {'^\Hqp\iI>e') , as in eq. (30 1 



(E~€)a^ dE'A{E')b{E'), {E - E')b{E') ^ A* {E')a. (41) 
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The existence of the homogeneous solution in the continuum at E' = E leads 
to the general result in the form of the principal value plus the delta-function, 



P.v. 



1 



iE)S{E-E') A*{E') 



(42) 



E-E' 

where the residue factor r(E) is real and has to be determined self-consistently, 



while in the analogous situation in eq. (25) we could simply use the standard 



identity (26 1. The first equation of (41) determines now 

-A{E) 



r{E) 



E 



\A{EW 



(43) 



where A(i?) is the simplified form of the principal value integral (29 1. Introduc 



ing the coordinate representation of the scattering problem, one can show [16] 
that this residue defines the phase shift 5, 



i&Yi6{E) 



r{Ey 



(44) 



Finally, the normalization of the full wave function by the delta- function of 
energies in the continuum, allows one to find the remaining amplitude 



\A{EW 



[£;-£- A(s)]2 + r2/4' 

of resonance shape with the energy shift A(i?) and the width 

T{E)^2n\A{E)\\ 



(45) 



(46) 



in agreement with the definition of eq. ( |30[ ). Both, the energy shift and the 
width, are energy-dependent. The shift A(_E) vanishes if one can neglect the 
energy variation of the amplitudes A{E). The physical characteristics changing 
rapidly as a function of energy near the resonance peak lead to the highly asym- 
metric curves for cross sections which is a well known feature of autoionization; 
the typical example of the 2s2p excited state in the helium atom was analyzed 
in detail in the original work |16j . 



3.4 Overlapping resonances and super-radiance 

Our next step is the analysis of the eigenvalues and the eigenfunctions of the 
effective non-Hermitian Hamiltonian T-L. To emphasize the physics of super- 
radiance and related transformations of the resonance spectrum, we neglect for 
a while the energy dependence of the amplitudes A^. This assumption will 
be lifted in the continuum shell model when the threshold dependence of the 
amplitudes will be important for the spectroscopy of the resonances. If our 
energy is not close to thresholds, the energy dependence is smooth. In the 
case of many resonances with the same quantum numbers, such as the low- 
energy neutron resonances in heavy nuclei, the energy dependence is common 
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for all of them and can be usually subtracted [57] using the so-called reduced 
widths. Special precaution should be taken if the compound resonances come 
from the single-particle resonance that served as a doorway state on the way 
to a compound nucleus In this case the shape of the doorway strength 

function may distort the reduced widths (we return to the width distribution in 
Sec. 7.2). 



The two parts of the effective Hamiltonian (27) do not commute in gen 



eral and cannot be brought simultaneously to a diagonal form by an orthogo- 
nal transformation of the basis. The complex diagonalization of T-L leads to a 
biorthogonal set {XnXr} of right and left eigenfunctions which have complex 
conjugate eigenvalues £r and £* . All eigenvalues are located in the lower part 
of the complex energy plane; the resonance widths r^, eq. ([6|, are positive. 
This can be easily seen from the properties of the effective Hamiltonian (27). 
Its eigenvalue equation reads 



QQ 



1 



QV 



VQ 



Taking the matrix element with the biorthogonal partner Xri obtain 



QV 



1 



1 



H 



vv 



vv 



H 



VQ 



X; 



(47) 



(48) 



The principal value terms cancel, and the imaginary part of the eigenenergy 
determines the positively defined width (the so-called Bell-Steinberger relation): 



Tr = 2Tr{x:\HQrdiE - Hrv)HvQ\Xr 



(49) 



It might be convenient to work instead with a complete orthogonal basis 
of one of the parts, H or W, which will be distinguished as the internal and 
doorway representation^ respectively. Let \a) be the eigenfunctions of H with 
real eigenvalues [as the Bloch waves in the case of the lattice, eqs. ([s]) and ([9|], 



H\a) = ea\a). 



(50) 



The part W of the total Hamiltonian remains factorized after the unitary trans- 
formation to the basis a), compare eq. 

W^p = 27rJ2 {c^miA^ (2|/3) ^ 27r ^ B^^B^^ . (51) 

l,2;c c 

In the internal representation, we deal with the Hamiltonian 



c(opcn) 



According to eq. ( 34 1 , the complex eigenvalues are the roots of 



Det <^ 1 



k{£) 



(52) 



(53) 
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Two opposite limiting regimes naturally come into consideration. At weak 
coupling to continuum, the anti-Hermitian part is a perturbation to a standard 
intrinsic Hamiltonian of the shell-model type. Then we expect the appearance 
of narrow resonances corresponding to long-lived states barely coupled to decay 
channels. At strong coupling, the main role is played by the part W . In general, 
the whole picture including the resonances depends on the running energy E of 
the system, although this might be not important far from thresholds. 

Let us make simple estimates assuming that the intrinsic states \a) are ap- 
proximately of the same degree of complexity and have a similar strength of 
continuum coupling so that it makes sense to use the average parameters. The 



transformation ( 51 ) preserves the trace of W which we denote w. This trace will 
be the same after the final complex transformation to the eigenbasis of T-l, and 
the typical width of a resonance in the case of weak coupling will be F ^ w/N. If 
this width is small compared to the spacings D between the levels , the usual 
perturbation theory should be valid. This can be expressed with the aid of the 
coupling parameter k ^ (T/D) <C 1, or, in global characteristics, w/{ND) ^ 1, 
the total summed width is small compared to the energy interval covered by the 



relevant intrinsic levels. Only the diagonal part of W in the representation (52 1 



is important in this limit, and the resonance states, as in the example (15), still 
can be labeled by the quantum numbers of H, eq. ( |50| , 

=e„- ^r„, T„ = 2nJ2\K\''- (54) 

c 

In this regime the reactions reveal narrow isolated resonances la) with partial 
widths 

7^ = 2^|i?^p, r„=^7S. (55) 

c 

This is, for example, the typical picture for low-energy neutron resonances. 

The situation drastically changes when we approach the strong coupling 
regime, k ~ 1. Then w ^ ND, the resonances start to overlap and the situation 
of similar isolated resonances turns out to be unstable with respect to the sharp 
redistribution of widths. To illustrate this point we start with the one-channel 



case, when eq. ( 53 ) reads 



a 

In the same basis we now can see a very different solution in the limit of k 3> 1, 
namely an emergence of a very broad state that accumulates almost the whole 
summed width w — J2a^a- If w ^ ND, all intrinsic states would be covered 
by the total width, so that introducing the centroid of real energies weighted 
with their original widths 7^, 



^, (57) 

w 
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we immediately find from eq. (56) tlie broad pole £i = Ei — {i/2)Ti with £i 



close to the centroid (57) and the width Fi close to w. The remaining — 1 



states are in this limit very narrow (trapped). With corrections of the second 



order in 1/k ^ 1, this special root of eq. (56) is given by 



E.^e + ^Y.^J-'^^-f^, (58) 



T,^w-Ay ^"(^"7'^ . (59) 

ct 

The trapped states have narrow widths Ft ^ w/{Nk'^) which are much smaller 
than the average starting width 7 w/N. The more detailed analytical so- 
lutions can be derived in a model |35j with equidistant intrinsic spectrum and 
weakly fluctuating 7^. 

This example shows a generic mechanism of the formation of the super- 
radiant state. When the widths of individual resonances overlap and k > 1, the 
dominating interaction between the resonances proceeds through the common 
channels in the continuum. Roughly speaking, energy uncertainty of an unstable 
state makes a virtual decay of this state possibly end with the return in the Q- 
space to a neighboring overlapped state. This coupling leads to the resonance 
collectivization and the segregation of a super-radiant combination, the story 
one can identify in the two examples discussed in Sec. 2. It is easy to see that 
the super-radiant state |$i) is a superposition of intrinsic states best adjusted to 
the decay channel. Expressed as a vector in Q-space, it is nearly aligned along 
the decay amplitude vector B = {-Bq}, its transverse components are small, 
^ 1/k. In the example of a periodic lattice of Sec. 2, this corresponded to the 
localization of the super-radiant states at the spatial edges of the chain opened 
for the decay. The trapped states return to the non-overlap regime and belong 
mainly to the subspace orthogonal to B. 

We can also establish the hierarchy of time intervals corresponding to various 
parts of the dynamics involved in this transition. At k S> 1, the longest time 
is the lifetime of the trapped states, Tt ~ Ti/Ft ^ K{h/D) 3> fi/D, where 
the Weisskopf time ^ h/D [7D] describes the recurrence of the intrinsic 
wave packet in its evolution in the closed space without decays. The inequality 
Tt > Tw means that the wave functions of trapped states can be considered fully 
equilibrated as corresponds to the proper compound stage. The next time h/ND 
is even shorter than the recurrence time representing essentially the lifetime of 
a certain intrinsic basis state |1) with respect to the intrinsic spreading. Finally, 
the lifetime of the super-radiant state is the shortest, ~ h/Ti ~ {h/NDn). 

Going to the case of two open channels, a and b, we will have two A^- 
dimensional amplitude vectors, B° and B**. In accordance with that, in the 
regime of strong continuum coupling we find two broad states, \^±), which 
divide among themselves the lion's share of the total width w. Their complex 
energies, £± = E± — {i/2)r± depend on the amplitudes, 7°^'' = 27r|B°'''p, of 
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the vectors B and the angle between them |T^. In particular, after the super- 
radiant transition, k 3> 1, 



If the decay amplitudes are close to orthogonality, the channels become effec- 
tively decoupled and we obtain two super-radiant states (the left and right edge 
states in the lattice example). In the "parallel" situation, — 0, the state 
l^"-) with destructive interference of the two decay channels joins the set of 
trapped states, while the constructive interference leaves only one super-radiant 
state |$+) with r+ w. The long-lived N — 2 states are concentrated in the 
subspace orthogonal to that spanned by the vectors B^^''. 

We can qualitatively anticipate the results for a multi-channel case. The 
matrix W has a rank M equal to a number of open channels. This defines the 
number of its non-zero eigenvalues; for M = 1, the only nonzero eigenvalue of 
this matrix equals its trace w. Therefore, after the transition to the strong con- 
tinuum coupling, the natural basis would be that of the doorway representation, 
where M eigenvectors of W with non-vanishing eigenvalues are used along with 
N — M vectors from the complementary subspace. Here we expect the appear- 
ance of M broad states sharing the total width, while N — M states are getting 
trapped. Such a picture was clearly observed in the numerical calculations [5] 
of nuclear structure of ^^O in the region of the dipole resonance studied as a 
function of the continuum coupling strength. However, if the multidimensional 
angles t^"^ between the amplitude vectors fluctuate, the possible destructive in- 
terference can suppress the manifestations of super-radiance, especially in the 
situation of many open channels. 

The doorway physics will be discussed in more detail later. Here we just 
mention one simple but instructive model |52], a set of n single- fermion levels 
Cj, including the upper one, j = 1, in the continuum with complex energy ei — 
{i/2)j. In the system of fc < n fermions, the interaction mixes the configurations 
partially populating the level ei and opening the decay for other states. If all 
n\/[k\{n — ky.] many-body states have energy higher than the decay threshold, 
in the presence of the interaction they acquire the widths and shift down from 
the real energy axis into the complex energy plane. However, with the increase 
of the width 7, complex plane trajectories of some states turn back and again 
approach the vicinity of the real axis while the rest of states have monotonously 
growing widths. Fig. 1. The number of the latter (super-radiant) is (n — 
1)1/ [{k— l)l{n — k)l], the number of configurations of remaining {k— 1) particles 
compatible with the occupied doorway state |1) (these particles share the orbits 



The segregation of super-radiant and trapped states occurs in fact rather 
sharply when the parameter k reaches the value close to 1. This is illustrated 
by Fig. 2 taken from Ref. [12] where the system of 1000 intrinsic states 
coupled uniformly to 250 channels was modeled. The random distribution of 
the widths at weak coupling evolves releasing the short-lived states until the 




(60) 



J + l). 
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cloud of super-radiant states sharply separates from the rest of the system that 
contains trapped states with prolonged lifetimes. 
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4 Nuclear continuum shell model 



As mentioned in the previous Section, the story of the effective non-Hermitian 
Hamihonian was essentially started with the search for an appropriate theory 
for description of resonance nuclear reactions. Although we will argue that 
the theory is much more universal and will show the examples of its broad 
applications, still nuclear physics provides probably the most explored arena for 
using this approach. At least partly, this can be explained not only by historical 
reasons. The main trends of the development of modern nuclear science clearly 
lead to the problems, where the continuum effects are exceedingly important. 

Experimental nuclear physics shifted its center of interest to the nuclei far 
from the valley of stability. The new physical problems combine here the ques- 
tions of limits of nuclear existence with the astrophysical mysteries of the origin 
of chemical evolution and supernova explosion. The standard shell model is con- 
structed to work in the discrete spectrum of nuclear, atomic or molecular states. 
The beta-unstable nuclei close to the drip lines are loosely bound, in many cases 
only the low-lying states (or just the ground state) are particle-stable. There- 
fore the conventional distinction between structure and reactions (the division 
between the subjects of the research as well as between the researchers) becomes 
obsolete and it is necessary to find ways to overcome this barrier. In fact, a sim- 
ilar situation exists for the entire class of open or marginally stable mesoscopic 
systems, including various micro- and nano-scale channels for transmission of 
quantum information. 

In loosely bound nuclei any real or virtual excitation is related to the con- 
tinuum part of the energy spectrum. The sharp discrimination between nuclear 
structure and nuclear reaction disappears and our theoretical description should 
cover both regions. In the classical problem of neutron resonances in heavy nu- 
clei we also are working in the continuum very close to the threshold of the 
discrete spectrum. But there (we return later to this class of problems) the 
details of the underlying structure are less important since the long-lived res- 
onances correspond to chaotic states of the compound nucleus, where one can 
use statistical approaches justified by random matrix theory. In nuclei close 
to the drip lines we are interested in details of weakly bound structure and in 
the spectroscopy of individual resonances. Therefore we need the generalization 
of the shell model that would allow the full understanding of this region and 
unified description of structure and reactions. 

Here the approach based on the effective non-Hermitian Hamiltonian be- 
comes extremely appropriate (along with other kindred methods, like for exam- 
ple the Gamow shell model [311 [33]). The signatures of super-radiance were seen 
(however not fully understood) in numerical simulations by Moldauer [71| . The 
first clear indication can be found in Ref . [8, , where the dipole excited states in 
""^^O were considered in a simplified shell model version, with decay probabilities 
(level widths) taken as variable parameters. The results sharply changed with 
the degree of openness of the system. Depending on the number AI of open 
channels, as the escape widths were increased, the formation of M very broad 
resonances was observed, with other states becoming more and more narrow. 



24 



The analogy to optical super-radiance and the explanation of these effects were 
suggested in Refs. [3 HO] ; the first results of applications to the nuclear shell 
model were given in the review article by Rotter [28]. Our discussion will be 
mainly based on the later development |7U UHl [SHI 131] ■ 



4.1 Two-level example 

We start with the illustration of dynamics on the borderline between discrete 
spectrum and continuum. We assume that the intrinsic, Q, space contains only 
two relevant states coupled to the single open decay channel through real ampli- 
tudes Ai^2- This simple example ,29. in fact has many non-trivial physical ap- 
plications [751 [731 IZH [111 [TBI IZZ] including experimental studies with microwave 
cavities [75]. It might serve [73] as the starting model for such a nucleus as 
3"'^Li8, where the open channel contains the residual nucleus ^Li and emitted 
neutron pair. For the pair we assume only one state ^Sq of relative motion. In 
some sense, this is an original Cooper effect ^80] where one pair can be bound 
on the background of other fermions even if the primary attractive interaction 
within the pair is insufficient for binding; the attraction in the vicinity of the 
nucleus is also helped by the exchange of the collective excitations of the core 

[mill]. 

The most energetically favorable wave function of the Cooper pair is a com- 
bination of available for the pair orbitals, for example Qpi/2 and lsi/2 in our 
simplified model (in reality, probably, it has more components). Both single- 
particle states are in the continuum since the adjacent isotope ^''Li is unbound. 
The effective Hamiltonian includes interactions of both types, the usual Hcrmi- 
tian V, such as pairing, and the anti-Hermitian, through the decay channel, and 



it can be written in a general form (27), 



\v- {i/2)A^A2 £2 - {i^Al ) ' ^^^^ 

where we assume that the real energies ei,2 are both positive (above threshold 
that is set at zero energy), and the factor \/2tt is included into the amplitudes 

^1,2- 

The complex eigenvalues of this Hamiltonian are formally given by 

ei + £2 - (V2)(7i + 72) 
2 

± ^^(ei - £2)2 + 4z;2 - llL±I^ _ ^[(e^ _ - 72) + ^vA^A^], (62) 

where 71.2 = A^2- This result is only formal if, as it happens in the realis- 
tic shell model, the amplitudes A12 are energy dependent and have to vanish 
at threshold. A typical complication in such problems is that the location of 
threshold is unknown a-priori. A consistent solution would require to consider 
simultaneously the daughter nucleus and maybe even the chain of decays as we 
will show later. 
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The behavior of complex energies predicted by eq. (62) is non-trivial 
173) . The Hermitian perturbation v leads to the standard level repulsion; in 
the absence of open channels, the lower root E- is pushed down crossing the 
zero level at = €162 and becoming bound. The intrinsic states are mixed 
and, if their widths were very different, after the mixing they become more 
uniform, the widths are attracted. Contrary to that, the imaginary perturbation 
in the second order makes the levels to attract each other, while the widths are 
repelled, and this is the beginning of the road to super-radiation. As this process 
is determined by the ratio of widths to the level spacing, it occurs immediately 
for degenerate levels, as in the ideal Dicke model. Indeed, for ei = £2 = e and 
V = 0, we obtain 

E+ = e, =e- ^(71+72), (63) 

the super-radiant state absorbs the whole width while the remaining state is 
stable although its energy is embedded into continuum. In the same way, for 
any number of degenerate intrinsic levels, the real Hamiltonian H is given by 
the unit matrix, and the energy splitting is fully defined by W that produces 
M non-zero eigenvalues. Another new effect is the possibility of crossing for 
the complex poles [83 , analogous to the well known behavior of two coupled 
classical damped oscillators [84 . 

The non-universal dependence of amplitudes Ai^2 [and even of real energies 
ei_2 if the principal value part A(_E), eq. ( |29[ ), is important] on running energy 
E makes the analysis much more involved [131 EH E] , since the complex poles 
move as a function of energy, and the position i? of a resonance on the real 
energy axis is defined by the solution of Re£{E) = E. The two levels in our 
problem are then found by diagonalizing different matrices. One should carefully 
take care of the vanishing width at threshold, so for example, when the lower 
level becomes bound, E- = r_ = 0. In the model of ^^Li, one can assume the 
usual energy dependence from the potential model for p- and s-orbitals. 



jp{E) = f3E^^^, 7,(S) = aVE. (64) 

In this way we obtain the self-consistent solution when, simultaneously with 
v'^ = epCsi we come to = r_ = 0. For special conspiracy of the parameters, 
the system also reveals a bound state embedded in the continuum. In general, 
bound states, resonances and energy dependence of the amplitudes are strongly 
correlated. 

The same effective Hamiltonian determines the resonance part of the scat- 
tering amplitude, 

= {E-£_){E-£+) ' 

where we correct a misprint in the corresponding equation of [29] . Only in 



the degenerate case of eq. (63), we see a clear Breit-Wigner resonance on the 
super-radiant state, 
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The dark state with F = is not seen in the scattering. Fig. 3 taken from 
P5] shows the typical cross sections for the positive, 3a, and negative, 3b, signs 
of the product A1A2 at certain values of parameters. The peak at low en- 
ergy in the case of the positive product and > £162, when the lower level 
is loosely bound, demonstrates a typical example of the threshold maximum of 
sub-resonant nature, whereas the remaining broad upper resonance at higher en- 
ergy is not shown. For the negative product A1A2, the presence of the threshold 
singularity is responsible for the finite value of the cross section at zero energy. 

4.2 Realistic calculations 

In all realistic applications, the truncation of both, V and Q, spaces is inevitable. 
With a truncated set of single-particle orbitals, the space of many-body states in 
the nuclear shell model, built, with full antisymmetrization required by Fermi- 
statistics, on these orbitals, can be naturally identified with our Q-space. The 
configuration interaction within this space is given by a set of matrix elements, 
traditionally of two-body but in modern versions also of three-body character 
[551 [55] . The solution process is usually based on a large-scale diagonalization 
of the Hermitian Hamiltonian matrix. The single-particle states can be conve- 
niently defined in a spherically symmetric basis being labeled by the spin-spatial 
quantum numbers m) and the isospin projection r. The diagonalization 

then leads to the many-body states with good quantum numbers of total spin 
and parity. The deformed mean field and corresponding rotational bands can 
emerge from the solution if the original space is sufficiently large. 

In this spherical shell model, the interaction matrix elements are defined for 
any pair of initial and final two- or three-body states with the same constants of 
motion — total spin quantum numbers (J, M), parity, and isospin if there are no 
forces breaking charge symmetry. The number of such allowed matrix elements 
can be large but still by orders of magnitude smaller than the dimension of the 
Q-space. The values of the matrix elements can be in principle theoretically 
derived from the vacuum nucleon and meson interactions fitted by few-body 
observables. In practice, these values have to be phenomenologically adjusted 
in order to reproduce thousands of experimental data. Here the theory has 
already a rich experience and many successes. Unfortunately, such an experience 
is essentially absent for the matrix elements in the H-pq+Hq-p class for coupling 
of the intrinsic space with the outside world. 

In the first stage of the development, usually a potential of the Woods- 
Saxon type and simple short-range residual forces were used in all parts of the 
effective Hamiltonian; the open channels were restricted to the one particle in 
the continuum [871 ES] • The single-particle resonances emerge from the internal 
shell-model bound states embedded in the open continuum. In a similar way, but 
with a more advanced model of intrinsic dynamics, it was possible to describe 
the inelastic electron scattering with the knock-out of a nucleon [89j . 

Since the effective Hamiltonian is energy-dependent, in general the position 
E,.{E) and the width Tr{E) of a resonance are parameters of the eigenvalues 
(Im of T-L{E) which depend on running energy of the experiment. The strictly 
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formulated theory is capable of predicting directly the observable cross section 
a°'^{E) for given entrance and exit channels. Then the convenient working 
definitions of the resonance parameters are the roots of obvious equations, E,. = 
Er{E = Er) and = TriE — E^). The shape of the resonance curve is evolving 
with the energy of the experiment being in general non-symmetric with respect 
to the position Er [57| and reveals a cusp required by unitarity at the threshold 
for opening a new channel [501 ED HI] ■ The cusps are well known experimentally 
both in nuclear physics and particle physics and clearly seen in exactly solvable 
models [93] and in full modern continuum shell model applications [SH [95] . 

The routine identification of the resonance width with the product of a 
spectroscopic factor (occupancy of a single-particle orbital) and penetrability 
is invalid being approximately applicable only in the regime of separated reso- 
nances. In a more general situation of strong coupling, the width results from 
repeated excursions to the continuum and back. This should be especially im- 
portant for many-particle processes as fission and fusion of nuclei [5H] . Partial 
widths for specific decay channels also can be calculated, and in the regime of 
overlapping resonances, the total width is smaller than the sum of partial 
widths [HI] while the equality holds for well isolated resonances. Another im- 
portant conclusion is that in general it is not correct to arbitrarily parameterize 
the set of resonances as poles of the S'-matrix. This can violate subtle unitarity 
requirements which lead to correlations of positions and residues at the poles, 
even for relatively well separated resonances TO: . 

In strongly bound nuclei, the spectroscopy of low-lying states can be de- 
cently described without taking into account the continuum effects, although 
the widths of resonances and the connection to the scattering problem are lost. 
In nuclei far from stability the most interesting physics is just on the border- 
line of the continuum and above. Here, the elements missing in the theoretical 
foundation of the shell model are to be added phenomenologically, using general 
physical arguments and available data. The situation is facilitated at low ener- 
gies and in light nuclei, where only few channels are open. Then a practically 
reasonable approach is to add a phenomenological potential description to the 
usual shell model inventory. In some versions of the shell model this is anyway 
required if the starting point was a harmonic oscillator field and proper radial 
characteristics are not defined by the model. 

4.2.1 Example: Oxygen isotopes 

A good example can be the shell-model description of the oxygen isotopes [3S1 
[5ni 131] • Of course there are promising developments of technically different 
approaches to similar problems of shell model accounting for the continuum 
effects [28l [97l [98] . Empirically, the last particle-stable oxygen isotope is ^"'O. In 
the trivial shell model, this would be a core ^^O plus eight neutrons completely 
filling in the levels 0^5/2 and The spin-orbit partner 0^3/2 is already in 

the continuum (the neutron width of the 3/2+ state in -'^^O is 96 keV), while 
the intermediate level lsi/2 is weakly bound, if the level scheme is built on the 
base of ^^O. At low energy, and due to the extra stability of the doubly magic 
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^^0, it is reasonable to limit our Q-space by this sd-shell model. The relevant 
decay channels at low energy are single-neutron and two-neutron decays moving 
the system down along the same oxygen isotope chain. As we mentioned, the 
knowledge of the decay thresholds is rather critical. To know them we have to 
study the spectroscopy of the daughter nucleus, and hence to look at the whole 
chain in its entirety, all neutron-rich oxygen isotopes from the root ^^O. The 
two-neutron channels can be restricted by the emission of the correlated pair in 
the ^Sq state. It is important to have in mind that for stable states, when the 
decay widths vanish, this approach automatically reduces to the routine shell 
model without continuum, maybe with renormalized parameters due to the shift 
operator A{E). 

In practice we need to add to the standard sd-model formulation, sec for 
example [99l 1100) new parameters describing the two neutron widths under an 
assumption that only the lowest states of allowed seniority (equal to zero or 
one) are populated in the daughter nucleus. For excitation energy above ~ 4 
MeV, one needs to include the final states with higher seniority. The widths 
for the emission of a d-neutron are supposed to grow from the corresponding 
threshold being proportional to the energy excess to the power 5/2 characteris- 
tic for i = 2. The coefficients are found from the neutron scattering problem off 
the ^^O nucleus. Even in the over-simplified framework where the residual in- 
teractions were limited to their main components, pairing and monopole terms 
[25], the results were quite satisfactory indicating the rich perspectives of the 
whole approach. We have to mention that in principle the effective shell-model 
interaction could be also complex in the presence of open decay channels; this 
problem was not solved yet. In calculations of reaction cross sections, the po- 
tential phases, especially Coulomb for charged projectiles, can be included along 
with the dynamics based on the coupling Hamiltonian H-pQ + Hq-p. 

The rediagonalization of the effective Hamiltonian for daughter nuclei ac- 
counts for the structural rearrangement of the remaining nucleus after neutron 
decay. Fig. 4 [30] shows the results of the continuum shell model for the chain 
of oxygen isotopes on the base of the shell model interaction |100| . There is a 
good qualitative and often quantitative agreement with available experimental 
data. Among interesting predictions, one can mention the ground state width of 
only 112 keV for the first particle-unbound even oxygen isotope, ^^O, where the 
Q- value for the two-neutron decay to the last bound isotope ^'^O is predicted to 
be only 21 keV. In this case only sequential two-neutron decays were considered; 
this should be extended for correlated pair decays. 



4.2.2 Some technicalities 

A favorable technical feature of working with the effective Hamiltonian is the 
absence of the need to use explicitly the complex scaling and related tools. One 
can still keep a diverse experience of the routine shell model and introduce new 
computational methods for the non-Hermitian Hamiltonian problem |94| . 

One of such technical improvements is the propagator calculation using the 
time representation of the propagator (33) for f — as a Fourier image for 
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positive time of the Hermitian evolution operator that, in turn, is expressed as 
an expansion over Chebyshev polynomials |101j . T„(cosa::) — cos{nx), 



oo 



H)"(2-<5„o)J„WT„(i/). (67) 



The main advantage of the representation is in the long-time dynamics where 
the Bessel functions J„(i) lead to the convergence important for the energy 
resolution in the following calculation of the cross sections, strength functions 
and level densities. Such an approach avoids direct diagonalization of large 
complex matrices. 



With the constructed Green function ( 33 1 at real energy and for energy- 



independent shell-model part Hqq of the Hamiltonian, one can solve for the 



full propagator ( [32| using the Dyson equation or its equivalent (34) (the so- 
called Woodbury equation) where the coupling with and through the contin- 
uum includes only the matrix in the channel space, usually of a much smaller 
dimension compared to the intrinsic shell-model space. The full Green func- 
tion Q allows one to study the damping of quasistationary intrinsic states with 
non-exponential decay of broad resonances. This approach, supplemented with 
addition of the collectivizing Hermitian interaction that is absent in the stan- 
dard sd-shell model, also gives a possibility to study giant resonances in loosely 
bound nuclei, a subject of great interest and contradictory predictions. In this 



approach the principal value term of self-energy ( 29 ) can be included without 
difficulty and, in practical applications [M], it improves the description of reso- 
nances seen |102| in neutron scattering from exotic oxygen isotopes. 
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5 Doorway states 



5.1 Dynamics based on doorway states 

In a complicated system, it is often the case that only a subset of intrinsic states 
\q) of space {Q} connects directly to the {V} space of channels. The rest of 
states in {Q} will connect to {V} only when they acquire admixtures of these 
selected states of the first type. The special states directly coupled to continuum 
are the doorways, \d), see for example |103) . They form the doorway subspace 
{T>} within {Q}, and the corresponding projection operator will be denoted here 
as v. The remaining states in {Q} will be denoted as \q) and their subspace as 



The full Hamiltonian in Eq. ( 19 1 can be decomposed now as 



{H-p-p + Hxip + Hpv) ■ (68) 



Note that the terms and H~ are missing because in accordance with the 
doorway hypothesis they are negligible. Also note that the diagonalization of 



the Hamiltonian part in the upper line of eq. (681 would provide the states \q) 
with the components of \d) mixed with the \q) states. The two last items in the 
above equation couple the doorway states, and therefore also the \q) states to 
the open channels. 

Doorway states play an important role emphasizing collective effects in de- 
scription of scattering and reactions with many-body systems. They are often 
serving as the connection between the continuum and the dense spectrum of 
quasi-bound states. As we will now see, the concept of the doorway is, in cer- 
tain situations, closely connected to the notion of super-radiance |41j . 

5.1.1 The case of a single doorway 

Let us consider the case when there is only one important doorway state \d). 



The matrix elements of the effective operator W, eq. (30 1, in the intrinsic space 
are now given by 

{q\W\q')=2TT {q\Hvv\c}{c\HT,v\q'} (69) 

c(opcn) 

with the doorway assumption: 

{q\Hvp\c) = {q\d){d\Hvp\c), (70) 
where {q\d) is the amplitude of the admixture of the doorway into the \q) state. 



Eq. ( 70 ) now becomes: 



{q\W\q') ^27T{q\d){d\q') ^ |(d|i?Dp|c) , (71) 

c(open) 
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where again we have separable matrix elements of the matrix W, this time 
irrespective of the number, M, of open channels. The doorway serves as a filter 
for the coupling of the {Q}-space to the space of open channels {V}. As a 
result, we find one broad state with a width 



q c c 



(72) 



This width is nothing but the decay width of the doorway into all available 
channels. 



5.1.2 When is this approach vaUd? 

As discussed above, the criterion of validity of the super-radiant mechanism 
is that the average spacing between the levels in {Q}-space is smaller than the 
decay width of a typical state "before" the super-radiant mechanism takes effect. 
This can be expressed, in the case of doorways, in the following way [31]. 

Consider the spreading width, F^, of the doorway state representing the frag- 
mentation of \d) into compound states \q}. If Nq is the number of compound 
states in the interval covered by the spreading width, their average energy spac- 
ing is 

(73) 

Before the super-radiant mechanism is turned on, the average decay width of a 
typical \q) state is 

rl = 2nY,\(q\HQr\c)\\ (74) 

C 

which can be estimated by 



(75) 
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These estimates result in 



^ « 4- (76) 

^1 ^ d 

We conclude that the requirement for the super-radiant doorway mechanism 
to be valid can be formulated as 

^ > 1, (77) 

the decay width of the doorway should be comparable or larger than its spread- 
ing width in order for the super-radiant mechanism to be valid. This condition 
is not always obeyed, but in some cases, as we will see, it is satisfied. 
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5.1.3 Isobar ic analog resonances 



An outstanding example, where the condition (77) is fulfilled and the super- 
radiant mechanism may take place, is given by the isobaric analog state (IAS) 
\IAS) |104j that is obtained from the parent state \tt) with certain isospin T 
by acting with the isospin lowering operator T_ that changes a neutron into a 
proton: 

\IAS) = const • r_|7r). (78) 

In medium and heavy nuclei, the IAS are surrounded by many compound 
states \q) of lower isospin T< = T — 1. The Coulomb interaction violates the 
isospin symmetry fragmenting the strength of the IAS over many states \q) and 
giving rise to the spreading width F^^g of the IAS. If located above thresholds, 
the IAS can also decay into several con tinuum channels with the decay width 



Tjj^g. In such cases the condition (77) is usually satisfied. For example, in 



^°*Pb the spreading width of the IAS is about 80 keV while the escape width 
is about 160 keV [104] . thus the ratio ^\as/^\as ~ 2- This situation is typical 
for other heavy and medium mass nuclei. In the Z ~ 40 region the above 
ratio is in some cases even larger than 2 |104| . The super-radiant mechanism is 
therefore relevant to this case providing a straightforward explanation why the 
IAS appears as a single resonance with the decay width of the doorway \IAS), 

{IAS\HQr\V) \ (79) 



5.2 Giant resonances 

Giant resonances in nuclei (or atomic clusters) present another example of sim- 
ilar physics |105l I106j . Usually, the giant resonances are discussed in terms of 
particle-hole (p-h) configurations with identical spin-parity quantum numbers. 
The coherent residual interactions form a correlated state that carries much 
of the transition strength of a corresponding multipole operator. The giant 
resonances are mostly located in the particle continuum decaying via particle 
emission to the ground and excited states in the daughter nucleus. 

Since the p—h giant resonance |G) is surrounded usually by a dense spectrum 
of 2p — 2 ft, and more complex configurations, the residual strong interaction will 
mix the state \G) with this background. Each of the resulting states denoted 
as \b) will contain the admixture (G|6) of the giant resonance. When above the 
threshold, the mixed states \b) couple to the continuum. If we assume that the 
dominant coupling is via the admixture of the giant state, then \G) serves as a 
doorway, and the matrix W is separable, given by 

{b\W\b') = {b\G){G\b') Y,{G\A\c){c\A\G), (80) 

c 

where (GIAjc) are the decay amplitudes for the giant resonance state. The 
matrix (6|VF|6') is again of rank one and the non-zero eigenvalue is equal to 

F^ = 27r^|(G|A|c) (81) 
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This is valid under the assumption that the energy spread of the background 
states \b) is small compared to their decay width, 



Tl = 2nJ2\{b\A\c) . (82) 



As before, this condition can be expressed as 

^ G 



>1, (83) 

^ G 



the spreading width of the giant resonance is smaller than its total decay width. 



When the condition (83) is not satisfied, so that the energy intervals between 



the background states are larger than their decay widths fJ, the situation of 
a single decay peak for the giant resonance might not hold. Still one could 
expect some bunching of background states into groups with spacing within the 
group smaller than the decay widths. Each such a group then can be treated 
separately using the super-radiant mechanism and appear as a single peak in 
the decay (or excitation) curve. One will then observe intermediate structure 
resonances in the giant-resonance energy domain. 

One should remark that for many giant resonances, such as the giant dipole 
resonance (GDR) or giant isoscalar quadrupole resonance, the above condition 



( 83 ) is not fulfilled |105| I106j , and the spreading width of the resonance exceeds 
the decay width. However, in continuum RPA calculations it was found that 
high lying giant resonances have substantial decay widths |105] . larger than 
their spreading widths. For example, the isovector monopole resonances calcu- 
lated in heavy nuclei had decay widths of the order of 8-10 MeV, larger than 
the spreading widths which are of the order of a few MeV. In these cases the 



condition (82) for super-radiance is satisfied. 

The GDR in nuclei not far from the valley of stability is shifted up in energy 
compared to the unperturbed excitation energy ~ Ihij oi p — h states. The 
GDR strength then saturates a significant part of the classical dipole sum rule 
although some strength still remains around the unperturbed energy. The situ- 
ation changes for nuclei with a large neutron excess farther away from stability. 
The fraction of dipole strength in p—h states that are not shifted increases due to 
the possibility of cooperative oscillation of the neutron excess against the core. 
If in an unstable nucleus some of these states are already in the continuum, they 
interact not only by the residual dipole-dipole interaction but through common 
decay channels. Then the situation reminds what we have discussed for the 
main branch of the GDR. As a result, the so-called pygmy-branch of the giant 
resonance |107[ I108j is formed being currently a subject of serious experimental 
interest. This soft collective dipole mode serves again as a doorway state. Such 
a mechanism of formation of the pygmy-resonance was predicted in [109] and 



considered more in detail in |110j . 

In the limit of small energy spread of the p ^ h states, the interplay and 
redistribution of dipole strength and decay width to a channel c between the 
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main GDR and the pygmy-resonance is determined by the "angle" between the 
multidimensional vectors A — {A'^} of decay amplitudes and d = {dqo} of dipole 
matrix elements connecting p — h excitations and the ground state. If these 
vectors are "parallel" , the soft mode loses its collective character, and the main 
GR peak concentrates both the strength and the continuum width. A situation 
of this type takes place in the case of the collective A-isobar— nucleon-hole states, 
discussed in the Section 5.2, when the pion matrix elements determine both 
the strength and the width [38l llllj of the doorway state. In the extreme 
"orthogonal" case the main peak would concentrate the multipole strength but 
the whole width would go to the soft mode so that the main peak would become 
a "zombie" state embedded in the continuum. The reality corresponds to an 
intermediate situation with the weight of the pygmy-branch rising closer to the 
drip line. In such cases it makes sense to speak about exit doorways [1121 1113] . 



5.3 Several channels and intermediate structure 

Let us turn to the situation when the \q) states are coupled to a number of open 
channels in the T'-space, a typical situation in nuclear physics. In addition to the 
elastic channel, the compound states in many cases can decay to excited states 
in the daughter nuclei. In this situation we have a number of open channels 
M that is still much smaller than the number of \q) states. In many cases it 
is not the same doorway that connects to all the open channels. In fact often 
the doorways are specified, and their number is equal to the number of open 
channels |41) . 

We should point out the difference between the situation mentioned here and 
the case of a single doorway for all the channels discussed above. In the latter 



case the matrix W was factorized, eq. (70), as Wqq' = aqU*, but the present 
case is more complicated because each channel has its own doorway, so that the 
matrix W is of the form 

M 

Wqq>^J2«'- (84) 

When the number \q) states is much larger than the number of channels Af , cor- 
relations are present in the matrix Wqqi . In this case the cross section will exhibit 
several resonances, with intermediate widths Ti^t such that Tg <^ Tint < Tg p , 
where Tg.p. is the decay width of a single-particle resonance. Such intermediate 
resonances were observed in nuclear reactions long ago [301 HT] and the cross 
sections were said to exhibit intermediate structure. 

In some situations, several doorway states \d) appear separated by energy 
intervals larger than their total widths. The intrinsic states \q) couple to these 
doorways locally giving rise to spreading of doorways. According to the doorway 
hypothesis, the \q) states couple to the continuum only via the admixture(s) of 
the local doorways. Then the matrix W separates into several disconnected 
quadratic blocks of dimensions x Ud, where is the number of states \q) 
contained within the spreading width of the doorway \d). For each such a 



sub-matrix, eq. ( 71 ) is valid, and the problem is reduced to that of the single 
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doorway. The resulting picture is a series of resonances separated in energy with 
decay widths equal to the decay widths of the individual doorways. 



5.4 Highly excited doorway state 

As excitation energy grows, many open decay channels are getting accessible. 
This region has its own specific features. A doorway state is singled out by a 
specific "signal" , either a quantum number different from that of the majority 
of background states, or by a character of intrinsic motion and, correspondingly, 
a type of external excitation. 

A nuclear example of the first kind is given by the IAS that is a state of 
isospin Ty = T that is embedded into a sea of states of isospin T< = T — 1. The 
Coulomb interaction and a small component of nuclear forces violating charge 
independence mix the IAS to the background. This spreading width of the IAS 
increases with excitation energy and then saturates |114| . If at this excitation 
energy the lifetime of the T> state with respect to the decay outside is shorter 
than the lifetime of the background r< states, the IAS preserves its individuality 
|115l I116j . One could expect that the purity of the IAS will rapidly be washed 
out at high density of the admixed states. As it turns out, just opposite is 
correct: after increased mixing, at some energy the purity of the IAS is restored 
as confirmed experimentally |117j . In fact, this phenomenon was predicted long 
ago by Morinaga [I18| and Wilkinson |119j . A qualitative explanation follows 
from plausible kinetic arguments [120] . The fully quantum- mechanical doorway 
theory [121] naturally justifies this conclusion. 

Let the doorway state \d) be immersed into a rich set of ^ 1 background 
states The state \d) with complex energy Sq = eo — (z/2)7o can decay into 
continuum (the corresponding width is 70) and to be mixed with background 
states by the coupling matrix elements which can be considered as real. 
In turn, the background states, after diagonalizing their intrinsic coupling, can 
at this energy decay into continuum through many, M ^ 1, complicated and 
essentially uncorrelated evaporation channels; their complex energies are e^. Be- 
cause of the complexity of those states and presence of many open channels, the 
coherent super-radiant collectivization of the background states is suppressed. 
The coupling of the doorway state with the background defines new quasista- 
tionary states \ j) , and, similarly to the usual theory of bright stationary states, 
the complex energies £j are the roots of 

£^-£,-J2^<^^0. (85) 

Each root j with energy Ej = Ej — {i/2)Tj carries a fraction fj of the doorway 
strength, fj = l. Substituting the evaporation widths by their typical value 



7cv and calculating the imaginary part of eq. (85), we find |121 



/, = (86) 

70 - 7ev 
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This result has a clear probabilistic, or kinetic, meaning indicating two routes 
of each fine structure resonance to the irreversible decay, 

=7o/,+7cv(l-/,). (87) 

The further analysis shows the resulting picture and its evolution with exci- 
tation energy from the viewpoint of reaction cross sections and time delays. The 
characteristic parameter here is the ratio 'jcv/D ^ Nj^v/^^ of the evaporation 
width to the spacing of compound states; the spreading width of the doorway 
state is ~ ND. Instead of super-radiance in this situation of many uncorre- 
lated open channels, when the background states overlap, the fluctuations in the 
excitation of the doorway state are washed away, see Fig. 5, and the depletion 
of the admixed background states of opposite isospin occurs faster than their 
population through the doorway mechanism. This justifies the predictions of 
isospin purification |1181 1119] observed in |117| . Such effects can be also inter- 
preted as related to the quantum Zeno effect |122| . In kinetic language of Ref. 
|120j , this stage corresponds to temperature of the compound nucleus exceeding 
a critical value when the total evaporation width of the admixed states becomes 
of the order of the spreading width of the IAS. In a similar formalism, analo- 
gous hierarchical effects were considered |123[ 1124] for chaotic quantum dots or 
atoms coupled to the electromagnetic modes of a cavity immersed in a bigger 
box. At random coupling, one can see mesoscopic fluctuations and oscillations 
of the survival probability. 

Another related effect is the so-called disappearance of collective strength of 
giant resonances excited in a high-energy heavy ion collision |125j . The analogs 
of giant resonances are expected to exist being built on all other nuclear states 
{Axel-Brink hypothesis) . However, beyond a critical excitation level, excited in 
heavy ion collisions fine structure states, according to the mechanism discussed 
above, decay in the continuum through incoherent evaporation channels before 
they are capable to transfer excitation to a collective degree of freedom: in a 
sense, the original excitation has no time to open the door through the doorway. 

5.5 Additional comments 

We have considered few examples justifying a generic qualitative picture of 
formation of doorway states and intermediate structures in various nuclear pro- 
cesses. It would be interesting to extend this approach to other mesoscopic 
systems (the application of the doorway concept to giant resonances in atomic 
clusters was made in |126| ). The notion of the doorway state has been applied to 
atoms and molecules [127) , to quantum dots |128) and fullerenes |129L I130j and 
recently even to a classical phenomenon of earthquakes [131] . In the physics 
of cold atoms |132j . the coupling of the internal states to the continuum via 
the Feshbach resonance influences profoundly the system affecting not only the 
imaginary part of the effective Hamiltonian but also the real part. This enables 
one to manipulate the properties of the system in a trap by changing the sign 
and/or the strength of the two-body interaction. The Feshbach resonance in 
this case serves as a doorway. 
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The traditional theory of intermediate structure in nuclear reactions |1331 
I134j considers usually statistical sequences in time of multi-step particle inter- 
actions with possible pre-compound (or pre-equilibrium) decay to the contin- 
uum. The present formulation emphasizes complementary aspects of under- 
lying physics, namely those that follow from the strict quantum-mechanical 
description of complicated many-body dynamics generated by the mean field 
and residual interactions. The approach based on the effective non-Hermitian 
Hamiltonian rigorously predicts the existence of different scales in the energy 
dependence of observed cross sections. These scales are defined by the intrinsic 
dynamics that combines collective (coherent) and chaotic (incoherent) features. 
The new collectivity responsible for the emergence of the doorways and related 
intermediate structure appears, under certain physical conditions, due to the 
factorized nature of corresponding terms in the effective Hamiltonian. This sep- 
arability is deeply related to the properties of unitarity of the scattering matrix. 
In this way we obtain a supplementary description of many-body dynamics. 
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6 Other applications 



6.1 Atomic and molecular physics 

Many-body physics of complex atoms has many features in common with nuclear 
many-body physics — shell structure, residual interactions, growth of complex- 
ity of excited states with excitation energy, very complicated — "compound" 
— wave functions of mixed states with several excited electrons, presence of 
the continuum etc. This structure and high level density are revealed, in par- 
ticular, in the photoabsorption spectrum, see for example, |135l 1136] , Typical 
autoionizing states appear already in the set of two-electron excitations. 

The whole super-radiant formalism was translated to the language of atomic 
physics in Ref . [35] . The physical difference from analogous nuclear phenomena 
is mainly due to the long-range character of the electrostatic interactions. In 
nuclear neutron resonances typical spacings between the resonances are greater 
than their widths, and the super-radiant collectivization requires higher energy 
above the neutron separation threshold, where the effect might be overshad- 
owed by the opening of new incoherent channels. In atoms the influence of 
threshold factors, ^ kR, is suppressed by the Coulomb effects, which facilitates 
the appearance of collectivization through continuum. According to [3S], the 
photoabsorption experiments will be mostly sensitive to narrow trapped states; 
in atoms the intrinsic collectivity similar to that generating giant resonances in 
nuclei is substantially weaker. But in electron scattering from ions one should 
observe broad doorway states and, with sufficient resolution, narrow trapped 
states as well. This difference is clearly seen in numerical calculations |35j . 

Experimental data mainly received with multiphoton laser excitation of 
atomic beams confirm the trend to the formation of trapped states with the 
autoionization rate suppressed compared to the predictions for a simple atomic 
configuration |1371 1138] . Ref. |139j provides a detailed review of this research 
area. 

The practical calculations of atomic spectra in principle can be done with the 
aid of the same shell-model machinery (currently called by the term "configura- 
tion interaction" coined in quantum chemistry) as in nuclei, especially because 
the interaction is essentially known. A new element comes from the presence 
in the same energy region, along with complicated compound states, also the 
Rydberg states of simple nature. The Rydberg spectrum is condensed near 
every ionization threshold. However, the large difference in structure makes 
the mixing between compound autoionizing states and Rydberg autoionizing 
states weak being in addition strongly suppressed by their very different radial 
behavior. The real spectrum consists therefore of two weakly mixing subsets. 

Interesting physics appears in the case of a simple (for example, hydrogen) 
atom in a strong uniform magnetic field |140) . The bound states in the con- 
tinuum can be found according to the results of our two-level problem. Sec. 
4.1. In the field around the critical value, 3^ = rn^e^c^fi^ ~ lO^T, the cy- 
clotron energy huj^i and therefore, the splitting of Landau levels, are of the 
same order as atomic Rydberg energy me^/fi^. The Rydberg series overlap. 
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and the sequences of Landau levels in the continuum interfere. This creates the 
situation of super-radiation and trapping, and, for some values of parameters 
(magnetic field), bound states in the continuum. It was stressed |141j that this 
phenomenon, earlier found in specifically adjusted potentials, is more general 
leading to segregation of fast decaying and trapped states. 

In atomic physics one can directly calculate the necessary matrix elements 
between the internal wave functions and continuum decay channels. An approx- 
imate statistical treatment [35] of the states J'^ = 4~ in the cerium atom showed 
that these matrix elements have the Gaussian-type distribution that leads to the 
Porter-Thomas distribution (PTD) of the resonance widths. Later, Section 7, 
we discuss the recent experimental indications and corresponding theoretical 
arguments for nuclear neutron resonances which show generic deviations from 
the PTD. In atoms the presence of two components of the spectrum may distort 
the picture. 

Radiationless and multiphoton transitions in atoms and molecules present an 
appropriate applications of super-radiant physics. The effect of super-radiance 
and trapping in molecular transitions in the approximation of degenerate in- 
trinsic levels coupled to the same decay channel was considered in 0B] and 
possible generalization to non-degenerate levels was mentioned (called in this 
paper "quasitrapping"). The actual presence of fast and slow dynamics in the 
multiphoton excitation was noticed long ago |142j . 

6.2 Intermediate and high energy physics 

As discussed in previous sections the super-radiant mechanism was applied to 
a number of different systems involving many-body physics. The phenomena 
discussed in the field of nuclear physics involved low energies. The universality 
of the super-radiant mechanism leads us to consider also the field of intermediate 
energy phenomena. The first discussion concerning this can be found in Ref. 
|37j . In this work the super-radiant mechanism was applied to a number of 
situations involving intermediate energy hadronic systems. The main objective 
of these studies was to point out that narrow resonances can be formed due to 
the super-radiant mechanism. Here we bring some examples. 

6.2.1 The NN states 

In current relativistic theories |143j of nuclear structure the Dirac vacuum of 
occupied negative energy states plays an important role in describing some of 
the nuclear properties. The influence of these negative energy states, N, is 
often implicit. It was suggested |144) that such theories could be tested by an 
observation of particle-hole states, in which a nucleon is raised from the negative 
energy Dirac sea and placed in an unoccupied state with positive energy. That 
we call an NN excitation. In most relativistic models the nucleon occupying 
a negative energy state experiences a strong potential of the order of several 
hundred, up to 700-800, MeV. Thus the excitation energy of the NN state is 
much lower than the mass of two nucleons, E^j^ < 2AI]y. The observation of 
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such sub-threshold nucleon-antinucleon excitations would provide a direct test 
of the relativistic theories of nuclei. A question should be raised about the 
feasibility to observe such excitations because of the expected large annihilation 
width. The NN pair is annihilated into mesons, in particular into several pions. 
One can expect that the NN states will couple strongly to the annihilation 
channel and decay very quickly. The large widths of these states will hamper 
their observation. 

It is worth noting however that conditions for the formation of a super- 
radiant state exist in this case |144j . The spectrum of negative energy states 
is within a given shell. (The spin-orbit splitting is nearly zero for the negative 
energy states). The unoccupied positive energy states form also a dense spec- 
trum, and therefore one expects a strong degeneracy among the unperturbed 
NN configurations. The energy splitting between various NN states is typi- 
cally 5-10 MeV, thus smaller than the individual decay widths into the common 
annihilation channel. We are therefore in the regime of overlapping resonances. 

Using a relativistic random phase approximation (RPA), a calculation of en- 
ergies and widths of NN states was performed in the framework of the cr — w 
model |145| for ^^O. The width was obtained from an imaginary term in the 
meson propagators that mimicked the coupling to the annihilation channel. The 
result of the diagonalization of the complex RPA matrices is quite amazing. All 
the states that resulted from the diagonalization had width of the order of 1-2 
MeV, while one state had a width of about 160 MeV. This is the super-radiant 
state. What happens here is that each NN state decaying to the common 
channel produces width of the order of 15 MeV. As a result of the diagonal- 
ization of the matrix, the wide super-radiant state emerges while the rest of 
states become very narrow and thus live a long time. One could anticipate the 
following situation. In a simple reaction, as for example nucleon pick-up, the 
super-radiant resonance, because of its large width will show up as a flat bump 
contributing to the background. On top of this background, narrow peaks will 
appear corresponding to the NN states that were "stripped" of their width by 
the super-radiant resonance. 

6.2.2 The A-Nucleus system 

Several years ago an experiment on inelastic electron scattering was performed 
with the reaction ^^C(e, e'p7r~)^^C. Evidence was found for narrow (several 
MeV wide) resonances in the A^^+^^C system in addition to a broad peak of 
approximately 100 MeV width jl46| . 

The free pion- nucleon resonance A33 with spin 3/2 and isospin 3/2 has a total 
width of Fa = 120 MeV, and the centroid of the four charge components is at 
1235 MeV; this is 300 MeV above the nucleon mass. The excitation of nuclei 
with various medium energy probes (electrons, protons, etc.) reveals a wide 
peak around the energy of 300 MeV. This peak is described in the framework 
of A particle- nucleon hole (AiV^^) configurations |111| . The A-component is 
treated as a particle bound in a potential that has similar characteristics to 
the usual nucleon potential. The nuclear A resonance is determined by the 



41 



coupling of the AA'^"^ configurations to tlie decay channel that involves pion 
emission. The A-nucleus system has properties that make the application of 
the super-radiant mechanism plausible. The AiV^^ configurations are strongly 
coupled to a single channel |c) — \A—1; {N + Tr)A), corresponding to a (nucleon 
plus pion)-system with quantum numbers of the A resonance and the residual 
{A — 1) nucleus. The partial decay width of each AA^~^ state is a fraction of 
the total width of the free A33, of the order of a few tens of MeV. On the other 
hand, the energy spacing between these states is of the order of several MeV. 
We are dealing here with the case of overlapping resonances. 

In Ref. [35] a simple shell-model description was used for the ^^C target 
nucleus. The nucleons occupy the lowest shell-model orbits Osi/2 and 0^3/2 
orbits. The A particle is assumed to move in a potential similar to the mean 
nuclear field. Therefore the energy spacings between the nuclear orbits and 
between the A orbits are close. One can then form AA^"-'^ states for various total 
spins J. For example, there are four A7V~^ states with J = 1+. Applying the 
theory of super-radiance we should find three very narrow states (with widths of 
several MeV) and one very broad resonance with a width around 100 MeV. The 
other states with spins 7 = 2+ and 3+ will add some narrow peaks. This picture 
resembles very much the observed cross-section [146j . One should bear in mind 
that other mechanisms might contribute to the observed widths of these quasi- 
stationary states. For example, the nucleon-hole states acquire spreading widths 
due to the residual nucleon-nucleon interaction. It is important to confirm the 
experimental results of Ref. |146| and try to study the possible narrow structure 
using alternative reactions that excite the A-nucleus resonances. 

6.2.3 High energy phenomena 

In many applications of the super-radiant mechanism, including the last two 
cases, the emphasis was placed not only on the wide super-radiant state, but 
also on the fact that the mechanism of the theory creates narrow resonances 
superimposed on a background formed by the broad super-radiant state, lead- 
ing in this way to the separation of fast (direct) processes from the long-lived 
structures. The consequences of this were examined and it was suggested that a 
mechanism of this type can explain the existence of narrow auto-ionizing states 
in atoms and narrow resonances in a number of strongly interacting hadronic 
systems. 

In the years 2003-2004 evidence was brought for the existence of a narrow 
baryon resonance with strangeness 5 = -t-1 around energy E — 1540 MeV 
|147l I148| . It was found in various experiments that the width of this resonance 
(coined as 0+ or Z"^), was a few MeV only. The resonance was situated on top 
of a very broad (over 100 MeV) background peak. It was difficult to understand 
the existence of this very narrow resonance structure in the framework of the 
usual mechanism for decay of baryonic resonances, whether exotic or non-exotic. 
The narrow width of the 5 = -t-1 exotic baryon indicates that the decay into 
the kaon-nucleon continuum is quenched either by selection rules, for example 
if the O"*" is an isotensor, or by some dynamical features. 
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In view of our previous discussions, it is possible that the narrowness of the 
peak observed is a result of the super-radiant mechanism "33] . This case serves as 
a direct example of the two-level super-radiant formalism discussed in Section 
4.1. In the simplest scenario one can consider the kaon-nucleon pentaquark 
system (su + ddu quarks) to be nonrelativistic and to form a quasi-molecule. 
Molecular-like structures for non-strange pentaquarks were discussed by lachello 
P_49J. The two particles interact via an attractive potential of a typical range of 
1 fm producing a p-wave resonance at energy 100 MeV above threshold. In Ref. 
|150j the authors estimate that the width of the resonance in such a potential 
is greater than 175 MeV, the value typical for strong decays of baryons in this 
mass region. A direct calculation gives the p-resonance width 190 MeV for the 
well of radius a = 1 f m obtained at the relevant energy; this required the depth 
of the well to be = 333 MeV; the width of 327 MeV is obtained for a = 2 
fm and 1/ = 14 MeV. In addition to quasi-molecular states, one should consider 
many-quark bag dynamics in a system of five quarks formed after the photon 
absorption by the original proton. Among the intrinsic states there are groups 
with the same quantum numbers, including "normal" quark states, quark-gluon 
states, paired states with singlet or triplet diquark(s), states with pions and so 
on. If the decay width of each individual state is greater than the energy spacing, 
the unperturbed spectrum will be that of overlapping resonances. Under such 
conditions, the dynamics will be completely dominated by the coupling to the 
decay channel, and the anti-Hermitian part of the Hamiltonian will be the crucial 
factor. As a result, the super-radiant mechanism will give rise to the observed 
spectrum with a very broad, ~ several hundred MeV, background peak of the 
super-radiant state and one or several very narrow resonances with width of 
order of few MeV. For example, for two overlapping resonances with AE' « 20 
MeV, and the bare width of 71,2 ~ 200 MeV, using the equations of Sec. 4.1, 
the estimate gives for the narrow resonance Fnarr ~ 4 MeV. A more detailed 
and quantitative account of this can be found in Ref. [39] . In Ref. |151) the 
authors proposed somewhat different candidates for the two interacting states. 
In summary, due to the overlap of unstable intrinsic states, the super-radiant 
mechanism may produce the narrow peak(s) on the broad background in exotic 
baryon systems with strangeness +1. 

We should remark here that the observation of pentaquarks is presently quite 
controversial. Some experiments that followed the initial claims of observation 
have reported the absence of pentaquark signals close to energy of 1540 MeV. 
Also some of the experiments that did report initially evidence for the existence 
of pentaquark have now retracted these claims [152] . Nevertheless the theoreti- 
cal discussion above is still valid. The assumption of the existence of two states 
close in energy around 1540 MeV might not be fulfilled. Generically, however it 
is possible to obtain a very narrow width for high energy resonances when there 
are favorable conditions for the super-radiant mechanism. 

In the last several years, in addition to the pentaquark, a plethora of pa- 
pers were published (mainly from the i?-factories) reporting the observation of 
new resonances jl53l 11541 1155j with small decay widths, much less than typi- 
cal hadrons would have at these energies (some of the reported resonances had 
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widths of a few MeV only) . Experimental and theoretical studies suggest that 
among these narrow resonances might be tetraquarks, that is systems with two 
quarks and two anti-quarks. These new resonances mostly lay around the mass 
region of 4 GeV, and the proposal is that they are composed of a light quark- 
antiquark pair and a heavy quark-antiquark pair, as for example ccuu, hhuu, 
and ccdd. Heavy tetraquarks are expected to decay into heavy states (as for 
example the J/1p^ charmonium) and light qq states, for example pious. The 
decay width of such processes is expected to be large, of the order of hundreds 
of MeV. However, in view of what we have discussed above, if there are two 
closely spaced tetraquark states with the same quantum numbers and both are 
coupled to the same final decay mode, there will be a rearrangement of widths, 
one state will become very wide while the other will become narrow. 

The most studied, so far, resonance among the ones suspected to be tetraquarks, 
is the X(3872) [HSl HSH HSg , which decays into (J/V')7r+7r-. The decay width 
is very small, F < 2.3 MeV. An unusually small width is an indication of 
strong incoherence and cancelation, suggesting that the super-radiant mech- 
anism might be responsible for this. The broad (super-radiant) state has a very 
large width and is embedded in the background. There are several observa- 
tions of additional resonances with widths not as small as that of X(3872), but 
still narrow compared to the typical characteristic decay widths in this energy 
region. 

The Belle collaboration reported recently the observation of the Z(4430) 
particle |156j . This particle is unusual because it has non-zero charge therefore 
providing a strong hint that it is a tetraquark. The dominant component of its 
wave function is probably ccud that decays into ^'tt. The decay width is of the 
order of 45 MeV. This again is an anomalously small value, especially considering 
the fact that the state is located several hundred MeV above threshold. It could 
possibly be viewed as a realization of the super-radiant mechanism. One should 
cautiously approach the new discoveries of tetraquark states and seek more 
experimental and theoretical confirmation before drawing certain conclusions 
about the mechanisms responsible for narrow resonance widths. The approach 
outlined here presents one of possible candidates to explain these long-lived 
resonances. 



44 



7 Quantum chaos in an open system 



7.1 Many-body quantum chaos 

Our discussion of doorway dynamics with separation of scales corresponding 
to external and internal interactions naturally leads to a general problem of 
quantum chaos. The deterministic chaos in classical mechanics is now a well 
understood and discussed in textbooks [157j phenomenon based on the extreme 
sensitivity of phase space trajectories to small variations of initial conditions. 
Contrary to that, quantum chaos |158j still is a subject of discussions, and 
some authors deny its existence recognizing only quantum signatures of clas- 
sical chaos |159j . In spite of this, we can formulate generic manifestations of 
quantum chaotic dynamics in a stationary quantum system as the presence 
of local correlations and fluctuations of observables approaching the extreme 
mathematical limit of canonical Gaussian ensembles of random matrices [160] . 

One-body quantum chaos emerges in systems of billiard type being caused 
by the boundary conditions, for example coming from an asymmetric shape that 
destroys spatial symmetries |161) . Such dynamics can be successfully reproduced 
experimentally with microwave cavities |162j , especially in the low-temperature 
superconducting regime with suppressed thermal dissipation |163) . In condensed 
matter systems, chaotic dynamics can be generated by the presence of disorder 
|164j or external fields |165j and, therefore, can be experimentally regulated. 

Our main interest is in physics of many-body quantum chaos where the 
chaoticity is originated by interparticle interactions [166] . As excitation energy 
and level density increase and any residual interaction mixing the independent 
mean-field particle configurations becomes effectively strong, many-body quan- 
tum dynamics in self-bound systems approaches chaotic limit. Then the sta- 
tionary states develop into exceedingly complicated superpositions of original 
simple configurations. The process of chaotization was tracked down in detail 
for complex atoms [167] . nuclei |168j and other mesoscopic objects. 

The ideal characteristics of quantum chaotic dynamics are given by the pre- 
dictions of the Gaussian orthogonal ensemble, GOE, (or Gaussian unitary, GUE, 
for the case of violated time-reversal invariance, and Gaussian symplectic, GSE, 
applicable, in particular, to time-reversal invariant systems with an odd number 
of fermions and therefore Kramers degeneracy) |1691 11601 11701 1171| . The distri- 
butions of energy eigenvalues in these ensembles are invariant with respect to 
corresponding (orthogonal, unitary, or symplectic) transformations of the basis 
and separated from the distributions of the components of eigenvectors, 

Pp{Ei, ...,En) = Cp{N) Yl \Em~Er,f exp[-/3NY,El/a 

m<n n 

Here Cfj{N) are normalization constants, /3 = 1,2,4 refer to GOE, GUE, and 
GSE, respectively, while a is the only available scale, a measure of the energy 
spread of the whole spectrum under consideration. Typical features of such dis- 
tributions are well known: level repulsion at small spacings (linear, quadratic 
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and quartic for those ensembles as can be understood from elementary quantum- 
mechanical perturbation theory [l72j ). small probability of large gaps in the 
spectrum, and the maximum of the probability close to the mean level spac- 
ing so that the levels in this limit form an aperiodic crystal with suppressed 
fluctuations. Counterintuitively, spectral fluctuations are much stronger for a 
regular motion when the level positions are uncorrelated and have the Poisson 
distribution. 

The nearest level spacing distribution that follows from eq. (87 1 is the first 
(and the weakest) feature of quantum chaos, in many-body systems it is reached 
already at a small interaction strength [168] . More sensitive measures deal with 
the fluctuations of the level number along large fragments of the spectrum [173] . 
Such an analysis is useful for checking the completeness and purity of long 
empirical level sequences, for example for low-energy neutron resonances which 
are important in practical applications |174l I175[ I176j . For the purpose of this 
review, we have to emphasize that such statistical applications assume that the 
long-lived resonances can be considered as representatives for the sequences of 
stationary states although in reality they belong to the continuum. 

The most interesting and far-reaching predictions of many-body quantum 
chaos refer to the structure of stationary states in the chaotic regime. Strongly 
mixed states in a narrow energy window have essentially the same degree of com- 
plexity |177j and the same macroscopic properties corresponding to thermody- 
namic equilibrium which provides a modern justification of statistical mechanics 
|166j . In the limit of random matrix theory, for large dimensions of the Hamil- 
tonian matrices, the components of the chaotic wave functions are uncorrelated 
Gaussian variables. Consider, for example, an experiment that singles out a cer- 
tain component of the wave function. Probing various states of similar structure, 
one would then discover that the distribution of corresponding probabilities is 
the chi-square for one degree of freedom, called the Porter-Thomas distribution 
(PTD) in nuclear spectroscopy |178j . The old experiments for neutron reso- 
nances confirmed this prediction. The joint analysis oi ly > 1 components, in 
a similar way, is expected to reveal the chi-square distribution for v degrees of 
freedom. 

The closer study of chaotic wave functions shows |179| that small perturba- 
tions acting in this region of states turn out to be enhanced compared to their 
effects for regular wave functions of simple configurations. This follows from 
the fact that, on the road to chaos, the matrix elements of simple operators be- 
tween chaotic functions of a similar degree of complexity decrease slower than 
the energy denominators in perturbation theory for the high level density. This 
chaotic enhancement is responsible for the unusually high, up to 10%, level of 
the parity violation effects in interactions of slow neutrons with heavy nuclei 
|180|, I181j . The chaotic mechanism of the enhancement is confirmed by the 
experiments on parity violation in nuclear fission by longitudinally polarized 
neutrons jl82| . The significant asymmetry of heavy fragments with respect to 
the neutron spin does not depend on the final distributions of product masses 
and kinetic energies being predetermined in the chaotic ("hot") stage of the 
compound nucleus. 
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Among other practically important applications of the ideas of quantum 
chaos we can mention the extraction of information about invisible fine struc- 
ture states in experiments with insufficient resolution |183j . New computational 
methods for large-scale diagonalization of Hamiltonian matrices of extremely 
high dimension are based on the chaotic properties of highly excited config- 
urations [1841 1185] and resulting exponential convergence of the quantities of 
spectroscopic interest — low-lying energies and other observables — as a func- 
tion of progressive truncation of the matrix preordered in a specific way. 



7.2 Statistics of resonances 

Many experimentally known sequences of complicated quantum states with iden- 
tical global symmetry in fact are resonances observed in the continuum spectrum 
and therefore having finite lifetime. A natural question is if the properties of 
quantum chaos predicted by the random matrix theory for closed systems are 
still valid for quasistationary states. 

Quantum theory of resonance states in complex systems is mostly devel- 
oped for two limits clearly distinguished, for example, in neutron reactions on 
heavy targets. Not very far from thresholds, the neutron resonances appear as 
quite narrow well isolated peaks in neutron cross sections. Each resonance here 
corresponds to a long-lived compound state probed in a specific mode of its 
excitation (entrance channel) and decay (exit channel). In elastic neutron scat- 
tering, the experiment singles out a certain component of the quasistationary 
wave function, namely the one corresponding to the continuum neutron (typi- 
cally in s-wave at low energies) and the ground state of the target. Here the 
resonance neutron widths, being proportional to the weight of this component in 
the compound state, are expected to obey the PTD, as for genuinely stationary 
states. In the opposite limit of broad resonances, their widths exceed the energy 
spacings. The cross sections in the region of overlapping resonances manifest 
a complex picture of wild interference {Ericson fluctuations |331I31]) with the 
correlational length of cross sections that is usually assumed to be determined 
by the average width of overlapping resonances. The latter quantity is defined 
by the inverse mean lifetime of complicated intrinsic states. The deviations from 
this limiting picture are known experimentally for a long time, see for example 
[Ml 11861 [55] and references therein. The extreme statistical averaging which is a 
core of the Ericson theory does not account for unitarity requirements which im- 
ply collective super-radiance and trapping phenomena occuring even if internal 
dynamics is chaotic. 

In the no-overlap regime of resonances in a chaotic many-body system, one 
can still expect the statistics of resonance parameters (energies and widths) to 
be in approximate agreement with the predictions of random matrix theory for- 
mulated for stationary states. This opinion was widely accepted but the recent 
accumulation of experimental data of higher precision [1871 1188j put this state- 



ment in doubt. Indeed, we mentioned, eq. (88), the stationary level repulsion 
at small distances. However, this is a property of a mixing by a Hermitian 
interaction. For unstable states, their energies are defined within uncertainty 
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related to their lifetime r ~ fi/T. There is no reason to expect the repulsion 
of resonances at short distances in energy comparable with their widths. The 
exact description in terms of the effective non-Hermitian Hamiltonian confirms 
the disappearance of such repulsion as can be seen already from our two-level 
model above, Sec. 4.1. Fig. 6 taken from Ref. |189j . where the case of a single 
open channel was studied, shows that, as a function of the increasing continuum 
coupling constant k, the probability -P(O) to have coinciding or very close po- 
sitions of resonances on a real energy axis (level crossing) is growing from zero 
and falls down again after the super-radiant transition when the broad state 
spreads into a part of the background. The dependence on the number of open 
channels and strength of the intrinsic interaction is illustrated by Fig. 7 |17] . 
In the GOE case all curves are qualitatively similar but the effect is getting 
more pronounced with the number of channels increasing. The model with two- 
body random interactions (TBRE), see also Section 7.4, shows that after the 
super-radiant transition P(0) falls down as a function of k in a universal way 
since, independently of the strength of the two-body interaction, we are dealing 
with uniformly trapped states. 

Moreover, the same qualitative argument of the preceding paragraph predicts 
also that energy centroids and widths of resonances cannot be independent, they 
have to be correlated, — the phenomenon clearly outside the standard theory. 
We need to construct a new statistical ensemble for unstable states generalizing 
the canonical Gaussian ensembles for closed systems. The similar ensemble of 
complex matrices [190) cannot serve as an appropriate model since its eigen- 
values are spread over the complex plane while the physical widths correspond 
to the negative imaginary parts of the poles of the scattering matrix — reso- 
nance widths Tj. are positive. This property is guaranteed by the construction 



of the effective Hamiltonian H [the choice of E'-'^'> in eq. (24)]. The simplest 
appropriate statistical ensemble for unstable states (our Q-space) assumes the 
intrinsic GOE dynamics [Hqq part) along with uncorrelated decay amplitudes 
(parameters of the coupling Hamiltonian H-po). These amplitudes include 
kinematic factors depending on energy and vanishing at thresholds. As men- 
tioned in our discussion of the continuum shell model, the energy dependence is 
crucial in realistic shell model studies. Dividing out this dependence we come to 
the reduced widths which can be considered as statistically independent. The 
PTD was expected to work just for these reduced quantities. 

Assuming the statistical independence and similarity of all decay amplitudes 
one can take them as random Gaussian quantities with zero mean and variance 

7" 



Aj^r = <5i2<5^^ ^, (89) 
where is a large dimension of the Q-space, and 7*^ parameterize character- 



istic partial widths. The ansatz (89 1 is natural for the system with chaotic 



internal dynamics. However, one can argue that, in realistic applications, the 



assumption ( 89 ) can be reasonable, independently of the presence of internal 
chaos, like in the case of weakly chaotic or even regular intrinsic dynamics, for 
N ^ 1. In self-bound rotationally invariant systems, such as a heavy nucleus. 
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we need to consider the states with the same quantum numbers J, M of angular 
momentum (plus parity and isospin). For a given shell-model partition, there 
are many ways of angular momentum coupling leading to the same total sym- 
metry. These paths expressed through the fractional parentage coefficients are 
essentially chaotic. The idea of random walk in the angular momentum space 
was used long ago by Bethe [191] for estimates of the spin-dependent nuclear 
level density. This geometric chaoticity plays an important role in statistics of 
self-bound systems |192|, 11841 11851 ^193] and it can justify the assumption (891 
even for the case of intrinsic interactions which are relatively weak to introduce 
the full dynamical chaos. Irrespectively of dynamic interaction, this leads to 
chaotic structure of wave functions and, for example, to the predominance of 
the scalar representation J = in the ground state of an even-even nucleus after 
averaging over all interactions allowed by the conservation laws [1941 1195| . 

With the GOE intrinsic dynamics and a single decay channel, it was derived 
algebraically [51 HH] that the common distribution function of complex energies 
([6| is given by the UUah distribution |196| . 

P{8„...,8^)^cJy[^] [n it'l l ^M-NF{E^,T^)], 

\ 71 ^ ^/ \m<n I ^ nl y 

(90) 

where the "free energy" of interacting resonances is given by 



^Er" + A Ern.r„. (9i) 



^ ^ — ' 7 ^ — ' 2a 



The distribution in the complex energy plane (90 1 reflects the interplay of 
various dynamical factors. The product of inverse square roots of the widths 
is the remnant of the PTD that is produced by the transformation from the 
amplitudes of given components in a complicated superposition of wave func- 
tions to their probabilities. The next product has the numerator responsible for 
the quadratic repulsion of energies in the complex plane, similar to the Gaus- 
sian unitary ensemble because the openness of the system effectively works as 
violation of time-reversal invariance. The denominator in this product shows 
"attraction" of the widths to their "electrostatic images" and therefore the ac- 
cumulation of the widths in the vicinity of the real axis (many long-lived states). 
The real axis remains a singular boundary that excludes the appearance of neg- 
ative widths. We can recall that the electrostatic analogy was fruitful in the 
first studies of quantum chaos by Dyson |197| . The free energy (91 1 is the 
main factor regulating the equilibrium arrangement of the complex energies. 
The super-radiance is contained in the last sum (repulsive width interaction). 
When typical widths grow beyond average spacings of real energies, the last 
term in F suppresses the probabilities of all configurations with uniform width 
distribution. The prominent maximum of probability is realized then on the 
configuration with one predominant width (super-radiant state). This means 
that the effect of super-radiance may coexist with chaotic internal dynamics. 

In many-channel situations, direct algebraic derivation of the joint distribu- 
tion function for complex energies is difhcult and the more advanced methods 
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related to supersymmetry should be applied. Such methods were developed 
first for closed systems |198) and then applied to unstable systems as well |199] . 
Unfortunately, the results are typically expressed in the form of cumbersome 
multi-dimensional integrals that makes the analysis not an easy task (the cal- 
culations are usually simpler for the intrinsic part described by the GUE rather 
than GOE). In this situation the direct numerical simulation of statistical en- 
sembles becomes the method of choice. We have seen in Fig. 2 [TT] the evolution 
of the distribution of complex energies as the strength of continuum coupling 
grows. The effect of super-radiance is clearly revealed as the segregation of 
super-radiant states (in the amount equal to the number of open channels) 
from the larger cloud of trapped states. This transformation is rather sharp 
and reminds a phase transition. Regrettably, the experimental data for multi- 
resonance systems with adjustable parameters that would allow to trace this 
restructuring are practically absent. The best candidates for this purpose are 
probably experiments with superconducting microwave cavities |163j and with 
acoustic chaos [200] . 

7.3 Statistics of resonance widths 

The statistical distribution of the resonance widths is an essential characteristic 
of a sufficiently long sequence of resonances. The low-energy neutron scattering 
off heavy nuclei provides rich statistics of such data. Usually the primary treat- 
ment of the data in the regime of well separated resonances is accomplished with 
the aid of the i?-matrix analysis [33] which we do not discuss here. At energies 
below the first inelastic threshold, neutron scattering can be only elastic, and 
we can apply the one-channel consideration. We will also omit a discussion of 
extracting the energy dependence that leads to the reduced widths. In practice 
this procedure should depend on the specific features of the compound system 
since the standard energy dependence defined by the orbital momentum of the 
neutron may be distorted by the shape of the single-particle strength function 
in the closed system [5^ . Assuming that we can have a reliable set of reduced 
resonance parameters, we can apply the statistical measures. 

For the set of physically equivalent reduced widths of a single-channel re- 
action one could expect the PTD. The admixture of intruder resonances cor- 
responding to different quantum numbers, such as p-wave resonances in the 
predominantly s-wave sequences, should violate the level spacing distribution 
and statistics of level fiuctuations. Such an impurity supposedly will move the 
distribution in the direction of the two-channel situation, and therefore increase 
the effective number of degrees of freedom given by the parameter v in the 
chi-square distribution, 

Puil) - Cj^/^-'e-^/^ (92) 

Here we can notice that even for well studied cases, such as ■^'^^U, the undetected 
admixture of p-wave can reach few percent |175) . Also we select an even-even 
target, so that all s-wave resonances have the same spin 1/2. For many years. 
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based on experimental data, it was thought that the PTD, eq. (92) with v = I, 
describes well the reduced neutron width distributions. However, recent precise 
measurements |1871 I188j claim that there are statistically reliable deviations 
from the PTD in the direction of the effective parameter smaller than one. 

The proper inclusion of the effects of coupling to continuum [201] with the 
use of the effective non-Hermitian Hamiltonian shows that such deviations in- 
evitably appear even for the sequence of statistically equivalent resonances al- 
though complementary effects of single-particle strength functions [69] can also 
be important. The PTD corresponds to the complete factorization of the joint 
distribution (87) into the independent GOE-like distribution of the resonance 
centroids En and the distribution of non-interacting widths. As seen from (861 
this occurs only in the extreme limit of vanishingly small widths r„ in compar- 
ison to the energy spacing D ^ a/N. The energy- width correlations enter the 
game and lead to the deviations from the PTD already at rather small value 
of the coupling strength k ^ j/D. The evolution of the width distribution 
as a function of increasing continuum coupling is illustrated by Fig. 8 [201] . 
Very soon, the chi-square model itself becomes a very bad approximation to 
the distribution predicted by (90). In agreement with the random matrix re- 
sults |202[ 1203] . the tail of the distributions behaves oc (P/P)"^ that cannot be 
described by chi-square distributions. 



7.4 Role of intrinsic dynamics 

Another interesting aspect of the width distribution problem is related to the 
intrinsic dynamics. The exact results cited above all refer to the canonical Gaus- 
sian ensembles of random matrices. In such cases of developed quantum chaos, 
all intrinsic states are of the same degree of complexity and their uncorrelated 
components give rise to Gaussian distributions, essentially as a manifestation 
of the central limit theorem. A slightly different type of statistics arises if the 
interparticle interactions, being still random, carry certain selection rules. We 
have already mentioned the role of geometric chaoticity in the process leading 
to complicated mixing of the nuclear wave functions even in the absence of 
interactions. 

Two-body interactions being the most relevant part of intrinsic dynamics 
bring a regular element in the Hamiltonian matrix of a closed system: many 
zeros and many repeated substructures |204j . Therefore, even if the resulting 
local pattern of the level repulsion and spacing distribution is in this case simi- 
lar to the GOE limit, the remaining correlations may have their influence upon 
the fine features, including the continuum coupling. This can be modeled with 
simulating the intrinsic Hamiltonian Hqq by the two-body random ensemble 
(TBRE). Here the two-body matrix elements are uncorrelated random quan- 
tities with zero mean and variance w^. The ratio v/d oi the TBRE strength 
to the appropriate mean-field level spacing characterizes the relative strength 
of intrinsic interactions. The estimates of the critical strength for the onset of 
chaos in the Fermi system of shell-model type with a certain number of particles 
in a truncated orbital space were given in Ref. |205| . 



51 



Fig. 8 shows the best value for the chi-square parameter v for different 
values of the continuum coupling strength k found |201j in the simulation of 
the dynamics by the GOE and TBRE with and without strong intrinsic interac- 
tion. The deviations from the chi-square distribution grow substantially faster 
for more regular dynamics both for one and two open channels. Intrinsic chaos 
is accompanied by the level repulsion and their relative ordering. Absence of 
chaos implies a stronger sensitivity to the continuum coupling. This physics sets 
the limits for the applicability of the traditional interpretation of the width of a 
neutron resonance as a measure for the strength of a pure neutron component 
in the wave function of compound nucleus. At the same time, we can acquire an 
instrument of evaluating the degree of intrinsic chaos by the observable width 
distributions for sufficiently clean resonance sequences. In general, various stud- 
ies allows one to conclude that increasing the intrinsic interaction we suppress 
the fiuctuations in the continuum making the compound states more uniform. 
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8 Signal transmission through a mesoscopic sys- 
tem 



8.1 Common physics of mesoscopic systems 

Scattering of an external beam off a many-body system is the most powerful 
instrument for obtaining the indispensable information on the internal structure 
of the system. At the same time, the situation can be inverted and our goal 
might be just the propagation of the signal itself |206| . The signal can be not 
only transmitted - we may need to transform it, split and recombine, separate 
resonance components, selectively amplify or suppress parts of the signal etc. 
It might be important to implement such different tasks simultaneously and 
record the results. All these tasks are necessary for quantum computing and in 
all modern applications of information treatment. A parallel problem is that of 
finding or synthesizing the quantum mesoscopic medium especially appropriate 
for a specific purpose. 

In fact, all these problems can be addressed in the framework of the effective 
non-Hermitian Hamiltonian. Here we have the intrinsic system with possible 
resonances along with open exit and entrance channels, plus some channels 
closed at a certain energy (frequency) of the signal but virtually influencing the 
results of a process. A system can also be embedded into complex environment 
and subject to external noises and decoherence but this set of problems is beyond 
the scope of the current review although it can be analyzed with similar tools. 

The terminology borrowed from nuclear reactions as a fully developed sub- 
field can be successfully used for other mesoscopic systems. In all cases of the 
signal transmission we are interested in the exit-entrance relation given by the 
scattering S'-matrix in the space of open channels. This matrix describes the 
probabilistic connection between the in- and out- quantum states. Without re- 



striction of generality, this unitary matrix is defined by eq. ( 38 ) , where we will 
be interested mainly in the second part that describes the signal propagation in- 
side the system (compound states) . The observable quantities are cross sections, 
or transmission probabilities, their correlations and Ructuations as functions of 
energy and intrinsic parameters of the system. 



8.2 Statistics of cross sections 

If the signal lives inside the system for sufficiently long time (longer than the 
Weisskopf time h/D, where 1/D is the typical density of states coupled to the 
signal), the situation becomes similar to the compound nucleus. In the experi- 
ment with good resolution we expect the cross section to reveal isolated narrow 
resonances corresponding to the levels of the compound system. As lifetime 
of intrinsic states shortens, the resonances broaden and start to overlap. Here 
one typically expects strong fluctuations of observable cross sections or trans- 
mission. The classical Ericson theory for the region of overlapping resonances 
|431 144j does not analyze the properties of individual resonances, concentrating 
instead on the observable statistical features. However, this theory does not ac- 
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count for the super-radiance effects since the couphng of intrinsic states through 
the continuum is not included into consideration. 

If we neglect direct reactions, the S'-matrix can be presented in an explicitly 
unitary form (38), 

S^l-^f='^mt, (93) 



where the ii'-matrix ( 35 ) is given by the Green function of intrinsic propagation 
sandwiched between the continuum amplitudes, 

K=A^^A^. (94) 

The scattering characteristics can be subdivided into average and fluctuational 
parts [207] . For example, it is assumed for the scattering amplitude that 

T'^^E) = T^^(i?) + T^^E), (95) 

where the averaging is performed over various resonances or, in practice, over 
an energy interval including sufficiently many resonances of supposedly similar 
nature. The same separation of contributions is performed for the ^'-matrix and 
cross sections. 

In the absence of direct reactions avoiding formation of the compound sys- 
tem, only elastic channels contribute to the average amplitude, T°-^{E) oc 5'^^ . 
In the language of the effective Hamiltonian, this is equivalent to the discussed 
above statistical ansatz (89 1 for the coupling amplitudes. Under this assump- 
tion, the averaging of the /-C-matrix gives ,46j 

sab a 



K-^{E)^~ib'''^K\ (96) 

where the density of states p ^ \/D should be taken in the center of the av- 
eraging energy interval containing N resonances. At ^ 1, this leads to the 
diagonal average S'-matrix, 

^ = (97) 

Assuming equivalent channels, k° — n, we can express the total cross section 
with the aid of the optical theorem, 



Ctot 



= ^a''^ = 2(l-ReS)=^, (98) 



where we omit the prefactor Air/k'^. 

It is convenient to work with the transmission coefRcient for each channel, 

-•^-l-l^-l'^ (T^T^- (99) 
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According to the definition ( 96 ) , this quantity, as well the average characteristics 



like (98), does not depend on the intrinsic dynamics if the parameter k is always 



renormalized for the proper level density. The so-called perfect coupling reg ime 



in channel a, r" = 1, corresponds to = 1 when the average 5'-matrix (961 



vanishes. This is precisely our region of the super-radiant transition. We can 



note that the transmission ( 98 1 is symmetric with respect to the substitution 
K 1/k; in the language of super-radiance this means that at large k one of 
N ^ 1 resonances broadens and disappears, and the transmission pattern in a 
given channel returns to the regime of no overlap. 

Similarly, the behavior of the fluctuating cross sections in various (by as- 
sumption, equivalent) channels can be encoded into another important param- 
eter, the elastic enhancement factor, the ratio 

F^^, b^a, (100) 

of elastic to inelastic cross sections |208j . In terms of the transmission coefficient 
and the elastic enhancement factor, we can derive 

Ft 



F + M -1 ^ ^ 

The conventional Hauser-Feshbach parametrization with cr°*' (x{l + S°-'') leads to 
F = 2. This result was repeatedly derived in various approaches j209| 11991 1210) 
in the regime of strong intrinsic chaos. The numerical analysis [3^ shows that 
the elastic enhancement factor goes down from _F w 3 to F = 2 with increase of 
the interaction strength reaching the standard value F = 2 in the GOE limit. 
This behavior requires more detailed understanding. Fig. 9 from Ref. [35] 
shows the average elastic fluctuational cross section in the TBRE simulation 
close to the critical coupling regime, n = 0.8, as a function of the intrinsic 
TBRE interaction strength A = v/d. We see the sharp transition to the GOE 
limit at the same value of A that corresponds to onset of chaos in the closed 
system. 



8.3 Fluctuations of cross sections 

It is accepted by the standard Ericson theory that in the overlapping regime 
the variance of the width distribution is small compared to the average width. 
Essentially this is based on the idea of partial widths being independent random 
variables. From our consideration of super-radiance it is clear that there is 
no room for Ericson fluctuations in the one-channel case because of a sharp 
transition from isolated resonances to the situation of one broad and many 
trapped states avoiding the overlapping regime. For many channels, Af ^ 1, 
one could expect that the partial widths are independent random variables, so 
that both, the variance Var(r) and the average total width F, grow linearly with 

2 

M. Then the ratio g = Var(r)/r falls down as 1/M. We have discussed earlier 
the width distribution for one open channel at realistic (small for the neutron 
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resonances) values of k. At k ^ 1, the ratio g is not small and strongly depends 
on K. As a function of M, Fig. 10, it falls much slower than 1/M revealing a 
trend to saturation; one can clearly see the dependence on intrinsic dynamics 
(parameter A) that is getting stronger as k grows. Again the onset of intrinsic 
chaos is smoothing the width fluctuations. 

In the case of intrinsic dynamics described by the GOE, the simulations 
agree with the prediction of a singularity at the critical coupling, 

Var(r) ln(l-K) 
.= ^oc^J--^. (102) 

This singularity is revealed, for any value of A, in the TBRE case as well. 
Another characteristic feature is the presence of correlations between the widths 
and the fluctuating part of the resonance amplitudes of the S'-matrix. Such 
correlations are amplified in the region of critical coupling. 

The conventional theory of overlapping resonances predicts that the corre- 
lation function of cross sections, 



C<,(e) = aiE)a{E + e) - a{E) , (103) 

as well as the similar correlation function Crie) for the scattering amplitude 
are essentially the same. Being expressed as ratios r(e) = C(e)/C(0), they have 
a simple form 

and the correlation lengths coincide, la- — It, being equal to the average width 
r. The qualitative argument in favor of this identification uses the idea that 
it is possible to keep correlations only on a time scale of the mean lifetime 
of the resonances, ~ h/T. Such arguments do not account for the effects of 
super-radiance. The analysis [IS] confirmed that the two correlation lengths 
are equal for large M, while 1^ < It for small number of channels, when the 
Ericson assumption are invalid, especially for weak intrinsic interaction. The 
fluctuations are not reproducible in different channels of the reaction. The 



correlation function is also different from eq. (104), in agreement with Ref. 
|211j . The equality Z = F is not valid, except for the case of very small k, as 
it was recognized long ago |186|. [571 1212j . Indeed, in the presence of the super- 
radiant trend, there is no justification for considering these quantities identical. 
Instead, the correlation length at large M agrees well with the expression [T^ 
through the transmission coefficient, 

I _ Mt _ M Ak 

Here the Weisskopf time Ti/D plays the role of the characteristic unit to measure 
the intrinsic life of the signal inside the system. For the same level density 1/D, 
the results do not depend on the intrinsic interaction strength A. The general 
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relation between the average resonance width and the transmission coefficient 
is given by the Moldauer-Simmonius formula |213l I214j , 



l--f^Hl-r). (106) 



Only at small r eqs. ( 105 ) and ( 106 ) lead to the equality I = T 



8.4 Prom Ericson to conductance fluctuations 

The so-called universal conductance fluctuations |2151I216T 11^1217] in electronic 
mesoscopic systems can be considered on the same footing as nuclear reactions; 
this analogy was repeatedly discussed in the literature [151 US |37] . The anal- 
ysis with the aid of the effective non-Hermitian Hamiltonian [l^ demonstrates 
how even delicate predictions of theory of universal conductance fluctuations 
are reproduced by the consideration based on a general quantum-mechanical 
foundation. 

The conductance of a mesoscopic system can be defined according to Lan- 
dauer as |42] a sum of partial cross sections, 

M/2 M 

G = E E (107) 

b=l o=M/2+l 

In order to establish the direct connection with conductance theory, we divide 
here our M channels, M considered to be an even number, into two equal groups 
of entrance, b, and exit, a, channels (this limitation certainly can be eliminated, 
see for example |218j ): we also omit the standard dimensional factor 2e^ /h. For 
equivalent channels, the average conductance can be directly expressed through 



the average cross section ( 100 ) and the elastic enhancement factor. 



G={ — \ a'^^^- — (108) 

V 2 y 4(F + Af - 1) ^ ' 

At M ^ 1, the limiting value Mr /A is fully determined by the number of 
channels and the transmission r. The simulations |48| are in good agreement 
with this prediction both for intrinsic interactions described by the GOE and 
by the TBRE. 

The fluctuations of conductance depend on the degree of chaos inside the sys- 
tem. The analytical results available for the GOE intrinsic dynamics |170| I206| 
|42] predict the variance Var(G) evolving between 2/15 in the ballistic regime 
to 1/8 in the diffusive regime. The approach of the effective Hamiltonian in- 
deed reproduces these values, see Fig. 11 taken from [15]. The crossover from 
one to another regime occurs precisely at the value of the intrinsic interaction 
parameter A that corresponds to the onset of chaos in the closed system decou- 
pled from the continuum. This and more detailed results on Var(G'), both for 
the GOE and for other models of internal dynamics, can be used for extracting 
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information on the degree of chaos inside the system from the experimental con- 
ductance data, for example, as a function of the degree of openness of quantum 
dots. Especially useful might be results on the various correlations of "elastic" 
(a = b) and "inelastic" (a ^ b) cross sections; the correlations turn out to be 
negative for the processes without a common channel. As a function of the con- 
tinuum coupling, the maximum sensitivity to the intrinsic chaos is observed at 
the intermediate strength of continuum coupling while at the perfect coupling 
the continuum effects overshadow the information coming from inside. An in- 
teresting object of application of these results can be found in correlations of 
speckle patterns PTOIE^ . 



8.5 Transport through nanostructures 

Without trying to give a full review of this extremely rapidly developing area, we 
only point out the relation of this tempestuous activity to the main topic of this 
article. In Section 2.2 we gave an example of the tight-binding model for one- 
dimensional lattice serving as a quantum wire for transmission of a signal, in this 
case propagation of an electron through a sequence of potential wells coupled 
by tunneling among themselves and through end barriers to the external world. 
This simplified model discussed originally in Refs. [iSj and [52j was recently 
generalized and put on a realistic foundation [IHl HO] ■ The whole approach with 
the idea of super-radiance and trapping in an open many-body system described 
by the effective non-Hermitian Hamiltonian was translated into the language of 
quantum dots in Ref. 221 . This is a prototype of many practical arrangements 
including semiconductor superlattices and quantum dot arrays |222l 12231 1224] , 
molecular bridges ^22?[|226] etc. Similar physics turns out to be appropriate for 
the description of energy transfer in biological "antennae" (like photosynthetic 
bacteria |227j ). The external barriers implement the controllable coupling to the 
continuum. It was shown [491 that a realistic sequence of potential barriers can 
be precisely modeled in the same approach and related to the open Anderson 
model [228| (the exact correspondence holds for weak coupling between the 
sites). 

A general situation with asymmetric coupling, recall the parameters and 
jfj in the Hamiltonian ([7]) and the solution for the resonances given by eqs. 



(12) and (131, was analyzed in detail in Ref. ^U\. Two realistically important 
cases are distinguished by the ordered character of the chain in contrast to the 
presence of internal disorder. In the case of a chain with N > 2 sites without 
disorder, the transport near the center of the energy band ([s]) reveals perfect 
transmission t = 1 at all resonance peaks for the critical value, 7cr — '^v^/q of 
the entrance width, where v is the tunneling amplitude, and we set 7^ = "f, "fR = 
"f/q introducing the asymmetry parameter q. The same value 7cr provides the 
maximum of the transmission integrated over the energy band. 



The transmission is determined by the sums ( 13 1 over the Bloch states inside 



the energy band. These sums can be calculated [50j as functions of e = E/2v 
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in the closed form, 



P±(e)^(l±l) 



v{z+ - z_) 



2N+2 _ -, 



(109) 



where z± 



1, a-i 



2N + 2, and a_ 



iV + 1. Inside the band 



e = cos/3 and z± = exp(±i/3). This leads to the "elastic" transmission (in this 
formulation back scattering) expressed in terms of rjLn = 7l._r,/2w, 



LL 



■ VLVR + Vr) + VR cos{2f3) + (1 + 77|j) sin /3 



('7L + '7i?,)(l + '7L??i?)(l + + 277L sin/3) 



(110) 



and analogously for 

by 



The " inelastic" channel (transmission) is characterized 



^VlVR 



(111) 



sin /3 

iVL + Vr){^ + VlVr) 

so that the condition of unitarity, r^^ + r^^ = -2lmT^^ is fulfiUed. The 
transmission integrated over the energy band of the width 4w and shown in 
Fig. 12 is well described by the same equation that is valid for two intrinsic 
states, 



(112) 



2v (g + l)[l+72/(4w2g)] 

The transmission maximum corresponds to 7 = 7cr, or ^lIr. — 4^2, the point 
of the "balance" between decay outside and tunneling inside the chain. This 
condition agrees with the results of Ref. |229j . 

The super-radiant effect, as we have already mentioned in similar context 
above, survives the intrinsic chaotic dynamics, in this case the diagonal site- 
disorder |52[ 150) . For weak disorder, the maximum transmission occurs at the 
same critical values of the continuum coupling. In the case of strong disorder, 
with the width of site energy fluctuation exceeding the tunneling strength w, 
the maximum transmission is shifted to higher values of the ratio 7/w. The 
width average reveals two super-radiant transitions at the expected values of 
the continuum coupling, ^ — 2v and 7 = 2qv. The details of this dynamics 
are discussed in ref. 50J, along with the comparison to the one-dimensional 
Anderson model [42l [230^ ,231j with symmetric and asymmetric leads. 

Coupled quantum dots have been already studied and used in the experi- 
ments directed to the development of quantum informatics [2321 1233] . The ef- 
fects similar to what have been discussed in Sec. 4.1, including the decoupling of 
special states from the leads that is equivalent to the appearance of bound states 
embedded in the continuum were discovered in |2341 1235j . There is a growing 
interest to the geometry more complicated than a simple one-dimensional chain 
|236l I237| . A three-dimensional transport diagram was probed |238j in triple 
quantum dots, where the conditions for stable transport were established. In a 
similar direction the search was performed with ^-junctions of superconducting 
Josephson chains 



|239j and T-shaped quantum dot structures [240) . Another 
possible geometry corresponds to a tetrahedral qubit |241j . 
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The theoretical results discussed above can be directly generalized to plane 
or three-dimensional networks of various structures |50| . It was already noticed 
p39] that the ^-junctions behave similarly to quantum Brownian motion on 
frustrated planar lattices. In the simplest models, continuum coupling is con- 
centrated at the perimeter. This arrangement can model an array of quantum 
dots or a particle tunneling through a lattice, with or without local disorder. At 
strong disorder, again the maximum transmission corresponds to the coupling 
strength growing proportional to the degree of disorder. Especially interesting 
can be the study of crossed chains that could be elements of future quantum 
computers |242j . The crossing point in such models plays a singular role con- 
necting all legs of the scheme. In addition to the Bloch energy bands, here, even 
with no disorder, the states of different nature appear, namely quasilocalized 
states with energy outside of the band and with wave functions exponentially 
decaying from the central point (evanescent waves). 

For the analysis of transport through realistic nanostructures it would be 
important to study, theoretically and experimentally, the influence of external 
noise (or "dephasing" |243j ) in addition to internal disorder. One inevitable 
source of noise is provided by electromagnetic fluctuations in the chain and in 
the background, see for example |244U245] . In our approach it can be modeled by 
adding many weak incoherent channels. Another extremely important features 
which should be introduced in the consideration are the effects of magnetic field. 
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9 Conclusion and outlook 



At the time of growing experimental and theoretical interest in open and marginally 
stable mesoscopic systems, it is useful to take a broad view of the whole field 
and find out the most general ideas and methods which would make possible 
to work out a more or less unified physical picture of various specific exam- 
ples. The approach related to the eff(X'tive non-Hc!rmitian Hamiltonian is quite 
general and at the same time easily adaptable to different cases. 

This approach, being in its starting formulation a rigoroiis application of well 
known quantum-mechanical ideas, keeps its generality and at the same time al- 
lows one to make proper approximations for concrete problems. The main new 
feature stressed by this approach is a collective interaction of internal states in 
a complex quantum system through their common continuum or radiation field. 
Being first suggested in quantum optics, this property turns out to be generic 
for all open quantum systems with many degrees of freedom. The process of 
such interactions can be traced all the way from almost independent well iso- 
lated resonances to super-radiance and trapping. Paradoxically, if the number 
of open channels does not increase too fast, the typical widths of narrow res- 
onances become even narrower with excitation energy increasing. We tried to 
demonstrate the unified physical picture of such phenomena in various subfields 
and such different objects as nanostructures and exotic nuclei, complex atoms 
and hypothetical long-lived heavy particles. The nuclear Ericson fluctuations 
turn out to be almost identical to universal conductance fluctuations in quan- 
tum dots; there are suggestions to look for similar physics even in such classical 
events as seismic waves. 

The unified formulation of the super-radiance physics of open systems al- 
lows physicists to understand better many questions which historically were 
asked and seemingly solved on different grounds and for different occasions. 
The segregation of broad and narrow resonances was first just an observation 
in complicated numerical simulations of nuclear reactions which seemed to be 
infinitely remote from coherent spontaneous radiation in the gas of two-level 
atoms. The details of Ericson fluctuations in the regime of overlapping reso- 
nances were not confirmed by the experiments — now we understand that the 
phenomenon of super-radiance (or perfect coupling to the leads in mesoscopic 
physics) related to the requirements of unitarity were not fully satisfied in the 
original formulation. We anticipate many such partial discoveries in further 
applications. 

As for any alive physical theory, we still have many open questions which 
hopefully will be addressed in the future. Our shortest list would include 

• clarification of the important problem of the role of internal chaos and 
interrelation between the regular or chaotic dynamics inside the system 

and its openness properties; 

• role of the multitude of weak external channels present in many practical 
cases, decoherence and equilibrium in open mesoscopic systems; 
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• better understanding of interplay between collectivity in intrinsic inter- 
actions, including superconducting correlations, and collectivity of super- 
radiant type through the continuum; 

• quantum signal transmission in complicated geometry; 

• development of more powerful computational methods. 
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Figure 1: A system of 4 fermions in 8 single-particle levels [8!/(4!4!) — 70 many- 
body states] with the upper orbital in the continuum. Particles interact via 
random Hermitian interaction. The lines (right panel) show the trajectories 
of many-body energies in the complex plane as the width of the upper single- 
particle level increases; 7!/(4!3!) many-body states become trapped and return 
to the vicinity of the real axis [52] . 
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75 

Figure 2: Segregation of the cloud of super-radiant states in the complex energy 
plane z = x + iy as coupling to continuum increases; here M/N — 0.25 [12]. 
Note the logarithmic scale on the vertical axis. 
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Figure 3: The near-threshold scattering cross section in the two-level model 
[55] 1 see text. The parameters of eq. (64) are selected from comparison with a 
corresponding potential model. 



76 



" r J III llHlttWfM-H- 




2Bf 



I 1 1 ■ 1 ' 1 ' 1 ' 1 ' 1 1 1 

■45 -40 -35 -30 -25 -20 -15 -10 
Energy (MeV] 



Figure 4: Continuous shell model calculations for the chain of oxygen isotopes 
[3D]. Bound states, resonances and cross sections of neutron scattering off the 
daughter nucleus are consistently calculated in the unified approach. 
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Figure 5: The relative probability of excitation of the doorway state from the 
continuum channel as a function of energy, E / D, and the evaporation width, 
7ev/-D, of the background states 112 1| for a specific value of the spreading width, 
r^/D = 10. The oscillations along the energy axis rapidly disappear as 7ev 
grows. 
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Figure 6: Nearest level spacing distribution found [189| by averaging 50 ran- 
dom matrices 160 x 160 with the Gaussian distribution of uncorrelated matrix 
elements and one open channel. The continuum coupling parameter k =0.1 
(panel a), 1 (panel b), 5 (panel c), and 50 (panel d). At critical coupling, panel 
b, the level repulsion at short distances disappears and the distribution can be 
fitted as a combination aPw{s) + (1 — a)PG{s) of the Wigner-Dyson distri- 
bution Pw{s) for the Gaussian orthogonal ensemble (GOE) and the Gaussian 
distribution Pg{s) centered at the mean value of s. After the separation of the 
super-radiant state, the distribution returns to the GOE regime. 
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Figure 7: The probability P(0) of finding the resonance energy spacing (in 
units of the mean spacing) s < 0.04 |47j. The intrinsic dynamics corresponds to 
the GOE, upper panel, with different numbers of open channels M, and to the 
TBRE, lower panel, for M = 10, and different strength of the two-body random 
interaction, from v = through the critical value for onset of intrinsic chaos to 
the strong interaction for degenerate single-particle levels, when the results are 
identical to those for the GOE case. Note that for the TBRE the vertical scale 
is different; at weak interaction the deviations from P(0) = start at a very 
small continuum coupling strength. 
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Figure 8: The best fitted parameter v for the chi-square distribution with v 
degrees of freedom used for the description of the width distribution found with 
averaging over 1000 random reahzations of intrinsic interactions for one, upper 
panel, and two, lower panel, open channels 07]. Full circles refer to the GOE 
intrinsic Hamiltonian, crosses stand for the TBRE with the interaction below 
onset of chaos, and squares for TBRE with no intrinsic chaos. The inserts show 
the unsatisfactory growing chi-square criterion of the fit as a function of k. 
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Figure 9: The fluctuating part of the elastic cross section for 10 open channels 
and K = 0.8 as a function of the intrinsic interaction strength in the TBRE; the 
GOE value is shown by a horizontal line HHI- The dashed vertical line shows 
the critical value Acr for the transition to chaos in the TBRE. 
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Figure 10: Normalized variance of the width as a function of the number of 
channels, for different coupling strengths k [15]. Circles refer to the GOE case, 
pluses to A = 0, crosses to A = 1/30, and squares to A — oo. While for weak 
coupling, K = 0.01, the variance decreases with the number of channels very fast 
in accordance with the expected chi-square distribution (dashed line) , for strong 
couplings. At = 0.5 and 0.9, the behavior is different from the 1/M dependence. 
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Figure 11: Variance of the conductance, eq. (107), for the TBRE model (7 
fermions in 14 orbitals with 20 open channels) as a function of the strength of 
intrinsic interaction 1481. 
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Figure 12: Integrated transmission, eq. (112), for N = 100 sites and different 
asymmetry parameters q [50] : squares refer to g = 1, circles to q — 10, and 
diamonds to g = 25. The vertical dashed lines indicate the critical value of the 
coupling strength in each case, where the integrated transmission is predicted 
to have its maximum. 
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